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CHAPTER  I 


INTRODUCTION 

Cercutiics  are  inorganic  non-metallic  brittle  materials  which  have 
been  used  by  mankind  for  centuries.  In  ancient  tines  they  were  used 
primarily  as  materials  for  pottery  and  artwork  because  of  their  unique 
properties  of  stability  and  aesthetic  appeal.  The  survival  of  pottery 
over  centuries  illustrates  one  of  the  great  advantages  which  ceramics 
possess  over  most  materials;  their  durability.  On  the  other  hand, 
ceramics  lacked  uniformity  and  reproducibility.  Until  three  decades 
or  so  ago,  users  of  ceramics  procurred  their  material  from  one  source 
and  one  particular  plant  of  a  supplier  in  order  to  maintain  uniformity. 
Ceramic  producers  were  reluctant  to  change  any  detail  of  their  pro¬ 
cessing  and  manufacturing.  The  reason  was  that  the  complex  material 
systems  being  used  were  not  sufficiently  known  at  that  time  to  allow 
the  effect  of  chcinges  to  be  predicted  or  understood. 

Today,  because  of  the  possibility  of  technological  control  of  the 
mechanical  properties  of  ceramics  emd  related  materials  they  have 
become  issues  of  major  importance.  Now  these  materials  are  used  in  a 
number  of  scientific  and  industrial  applications  such  as  turbine  blades, 
prosthetics,  ceramic  bearings,  and  the  like  [l].  The  rapid  growth  in 
particular,  of  the  electronics  industry  has  demanded  the  development 
of  many  specialized  types  of  ceramic  insulators.  One  of  the  important 
reasons  why  ceramucs  find  such  wide  application  in  a  wide  variety  of 
fields  is  that  these  materials  exhibit  a  wide  range  of  desirable 
properties.  For  example,  properties  of  individual  ceramics  include 
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resistance  to  heat,  oxidation,  corrosion,  and  abrasion;  high  elastic 
modulus,  high  strength,  and  low  density  with  applications  to  the 
nuclear,  optical,  magnetic,  electronic,  and  other  industries.  Other 
applications  are  much  too  long  to  list  here  in  detail,  but  some 
examples  will  serve  to  Illustrate  the  great  diversity  and  Importance 
of  these  materials.  Special  electrical  properties  form  the  basis  for 
applications  as  Insulators,  piezoelectric  transducers,  and  an  ever- 
increasing  family  of  semiconductor  devices.  In  the  past  twenty  years 
a  great  many  new  types  of  ceramic  Insulators  have  been  produced 
which  included  alumina,  magnesia,  beryllla,  zirconla  and  so  on.  The 
recent  use  of  barium  tltanate,  when  placed  between  the  plates  of  an 
electrical  capacitor  allowed  It  to  store  several  thousand  times  the 
amount  of  charge  that  could  be  stored  by  the  same  capacitor  if  air 
were  used  instead.  This  had  lead  to  the  production  of  a  series  of 
related  materials  called  ferroelectrics  which  have  helped  make 
electronic  equipment  lighter  and  smaller.  Magnetic  ceramics  Include 
cores  for  con^uter  memories  and  permanent  magnets  for  electric  motors. 

Special  chemical  and  thermal  properties  are  exploited  in  furnaces 
and  crucibles  which  are  essential  to  the  refining,  alloying,  and  cast¬ 
ing  of  metals.  ’Thermistor*  is  a  special  kind  of  ceramic  whose 
electrical  resistivity  varies  with  changing  temperature  and  is  used  as 
a  temperture  sensing  and  control  element  in  automated  manufacturing 
processes.  Thermoelectric  properties  are  important  in  developing  solar 
cells  supplying  energy  for  satellites  [2] .  The  high  elastic  moduli  and 
hardness  of  certain  ceramics  make  them  suitable  for  special  applications 
requiring  mechanical  stablity  such  as  gyroscope  mounts  and  wear 


surfaces  while  the  extremely  low  thermal  expansion  achievable  is  appro¬ 
priate  to  telescope  mirrors  and  to  applications  involving  thermal  shock 
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The  advent  of  the  missile  and  the  space  age  has  produced  exceptional 
demands  for  materials  which  will  stand  extremes  in  teiq)erature  and  in 
other  environmental  conditions.  Ceramic  materials  such  as  carbides  and 
graphites  are  being  used  in  some  rocket  nozzles  where  temperatures  in 
excess  of  5000°F  are  developed.  The  high  speed  of  sound  in  certain 
ceramics  has  led  to  their  successful  use  in  armor,  and  low  damping 
characteristic  of  many  ceramics  makes  them  useful  as  delay  lines. 

In  our  atomic  age,  as  atomic  energy  becomes  commercially  con^etltlve 
with  other  sources  of  energy,  more  efficient  nuclear  reactors  are  being 
produced  requiring  extreme  operating  temperatures  and  absolute  mechanical 
integrity  for  thousands  of  hours.  The  thermal  stability  of  certain 
nuclear  ceramics  has  led  to  their  adoption  as  fuel  elements  for  power 
reactors  and  in  many  cases,  several  functions  must  be  present  simultaneously. 
Because  of  their  extreme  hardness,  especially  at  high  temperatures,  ceramics 
are  widely  used  for  cutting  and  grinding  tools  in  Industry  [2] .  However, 
their  use  as  an  abrasive  has  been  limited  by  the  fact  that  they  generally 
do  not  form  sharp  cutting  edges  as  they  wear  away  (friability) .  Many 
different  optical  properties  of  ceramic  products  are  of  concern  in 
different  applications.  Perhaps  the  most  Important  are  those  optical 
glasses  and  crystals  used  as  windows,  lenses,  prisms,  filters  or  in  other 
ways  requiring  useful  optical  properties  as  the  primary  function  of  the 
material. 

In  citing  these  burgeoning  uses,  however,  one  must  face  the  historic 
problem  of  brittleness  and  unpredictability  in  mechanical  properties  of 


4 


1 


ceramics.  These  characteristics  still  remain  dictating  the  need  to  develop 
more  sophisticated  mechanisms  to  promote  general  design  with  ceramics. 

Further  a  polycrystalllne  ceramic  Is  usually  made  up  of  one  or  more  phases, 
grains  and  pores  of  various  sizes  and  distributions  and  Impurities  Inside  the 
grains  and  in  the  grain  boundaries.  All  of  these.  In  addition  to  process¬ 
ing  methods,  Influence  the  properties  of  ceramics  [3].  For  these  reasons, 
it  is  not  a  straightforward  task  to  formulate  design  procedures. 

Another  obstacle  to  the  wider  utilization  of  ceramics  Is  that  they 
fall  with  'glass-like'  brittle  fracture.  They  do  not  normally  exhibit 
appreciable  plastic  deformation  and  their  impact  resistance  Is  low  {4] . 

This  Is  one  of  the  important  reasons  that  the  application  of  ceramics  In 
an  engineering  sense  is  limited  by  these  relatively  poor  mechanical 
properties,  and  presently  the  only  generally  accepted  way  of  coping  with 
this  problem  Is  by  gross  overdesign  and  limiting  ceramics  to  areas  where 
structural  functions  have  secondary  consideration  [5] . 

If  ones  attention  is  focused  on  the  simplest  case  of  the  failure  of 
ceramics,  then  the  effects  of  multiaxial  stresses,  temperature,  strain 
rate,  and  time  dependency  are  neglected;  and  only  a,  the  uniaxial  stress, 
is  considered.  A  group  of  properties  —  fracture  toughness,  effective 
surface  energy  and  work  of  fracture  —  Is  also  important.  In  some  cases 
these  are  related  directly  to  strength.  This  group  gives  an  Indication  of 
whether  a  flaw  will  propagate  or  not.  Defining  Yi  is  an  effective  surface 
energy  for  the  initiation  of  fracture  from  an  inherent  flaw  of  size  c, 
a  statement  of  fracture  toughness  is: 


a 


(1-1) 
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where  Y  Is  a  geometrical  constant  and  E,  Young's  modulus.  (2Ey^)^  Is 
equal  to  the  stress  Intensity  factor.  Also  related  to  Is  the  work 
of  fracture  Yf*  defined  as  the  energy  required  to  generate  a  unit  area 
of  fracture  face.  In  general  Yf  ^ Yi  but  when  Yf  is  measured  In  a  test 
Involving  slow  controlled  crack  growth  (rather  than  a  partly  catastrophic 
crack  growth)  Yf  =  Yf  [6). 

From  the  engineering  point  of  view,  ductility  is  a  valuable  material 
property  not  only  because  ductile  materials  can  undergo  a  large  plastic 
deformation  long  before  failure  but  also  because  ductility  leads  to 
Increased  contact  area  where  high  stresses  may  occur  near  boundaries, 
dissipating  stress  concentrations.  Further,  ductile  materials  such  as 
metal  can  be  made  and  purchased  in  standard  forms  such  as  sheet,  tube, 
and  rod  etc.,  which  can  be  deformed  into  the  required  shape  and  joined 
to  form  structures.  Brittle  materials  must  in  general  be  used  in  the 
shape  in  which  they  emerg-.:  from  the  factory.  In  a  more  specific 
engineering  sense,  brittle  materials  have  the  characteristic  that  they 
are  unable  to  disperse  high  local  stresses  which  results  in  undesirable 
consequences  for  design  engineers.  The  very  low  ductility  and  strain 
to  failure  accompanying  brittle  materials  result  in  a  fundamental 
distinction  between  the  mechanical  strength  characteristics  of  ductile 
and  brittle  materials.  Material  property  test  results  for  nominally 
identical  steel  specimens  such  as  yield  and  ultimate  strength  will 
seldom  differ  from  specimen  to  specimen  by  more  than  5%;  and  a  minimum 
value  or  the  mean  value  alone  is  a  satisfactory  basis  for  design  whereas 
a  particular  'static'  strength  property  test  on  brittle  specimens  may 
show  a  variation  of  1007.  or  more  [7],  depending  on  the  number  and  size 
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of  the  specimens.  It  has  been  reasonably  established  that  this  variability 
is  not  a  result  of  badly  controlled  specimen  preparation  or  test  pro¬ 
cedure  but  it  is  an  inhereut  characteristic  of  brittle  materials  which  do 
not  exhibit  plastic  deformation  under  stress  [8].  Thus,  the  average  or 
mean  strength  value  is  not  an  adequate  basis  for  specification  for  engineer¬ 
ing  design;  the  variability  must  be  assessed  and  taken  into  account  in 
brittle  materials.  Designing  with  brittle  materials  does  not  require 
simply  different  numerical  values  for  the  properties  of  material  but  a 
new  concept  is  required  to  define  the  strength,  and  more  important,  the 
designer  must  think  in  terms  of  probability  rather  than  virtual  certainty. 

One  such  concept  is  known  as  the  ' weakest-llnk'  concept  and  it 
assumes  that  fracture  of  the  bulk  specimen  is  determined  by  the  local 
strength  of  its  weakest  volume  element.  The  'weakest-link'  hypothesis 
states  that  brittle  materials  fail  when  the  stress  intensity  at  any 
one  flaw  in  the  material  reaches  the  critical  value  for  crack  propa¬ 
gation.  If  flaws  are  uniformly  distributed  throughout  the  material 
then  the  number  of  flaws  present  in  a  specimen  is  proportional  to  its 
volume.  This  is  one  reason  why  large  specimens  are  expected  to  be 
weaker  than  smaller  ones  [9].  In  addition  to  these  variations,  the 
use  of  brittle  materials  in  structural  components  introduces  two  more 
difficulties  to  efficient,  reliable  structural  design.  They  are  (1) 
variability  in  mechanical  properties,  including  the  infrequent  but 
real  possibility  of  low  values,  particularly  in  large  volumes  of 
material,  and  (2)  the  accurate  definition  of  the  important  lower  proba- 
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bility  of  fracture  part  of  the  probability  distribution  curve,  which 
is  necessary  if  a  reliable  prediction  of  failure  probability  is  to  be 
made  for  large  volvunes  [lO] .  Since  catastrophic  failure  is  commonly 
induced  in  the  brittle  materials,  once  the  'fracture*  stress  has  been 
reached,  additional  and  more  information  may  be  needed  for  the  use  of 
such  materials  in  an  engineering  sense. 

The  weakest-link  concept  has  been  widely  used  in  developing 
various  statistical  theories  of  strength  which  differ  from  each  other 
only  in  the  way  in  which  the  use  of  the  concept  is  justified,  or  in 
the  assumed  form  of  the  distribution  fimction  of  local  strength.  Its 
application  to  a  solid  volume  was  first  proposed  by  Weibull  [ll]  who, 
arrived  at  a  distribution  function  which,  now  is  being  widely  used  for 
the  statisticcil  interpretation  of  a  variety  of  test  data  [l2j  .  This 
theory,  applied  to  brittle  materials,  addresses  the  statistical  varia¬ 
tion  of  strength  and  the  size  effect  due  to  flaws.  By  assuming  that 
only  tension  contributes  to  fracture  (fibers  in  compression  have  a  zero 
probability  of  fracture)  ,  Weibull  extended  the  theory  to  bending,  and 
to  any  uniaxial  combination  of  bending  and  tension.  Weibull  also 
extended  his  theory  to  failure  for  polyaxial  stress  states,  although  it 
has  been  questioned  by  many  authors  [13,14]. 

The  effect  of  a  biaxial  stress  field  on  the  fracture  behavior  of 
polycrystalline  alumina  ceramics  was  studied  by  Broutman  and  Cornish 
I is] .  The  various  stress  states  of  tension-compression,  tension-tension 


and  compress ion- conpress ion  were  generated  in  the  walls  of  thin-walled 
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cylinders  by  the  use  of  combinations  of  internal  pressure,  external 
pressure,  and  axial  loading.  Their  results  for  the  two  stress  ratios 
investigated  in  the  tens ion- tens ion  region  indicate  that  the  biaxial 
tensile  strength  is  less  than  the  uniaxial.  The  results  in  the  tension- 
compression  guadramt  were  obtained  by  a  combination  of  axial  compressive 
loading  and  internal  pressure  to  produce  hoop  tension.  The  experimental 
results  Indicated  that  the  tensile  strength  increased  when  compressive 
stress  existed  in  the  direction  normal  to  the  tensile  stress.  Weibull's 
theory  did  not  fit  the  data  well  as  large  amounts  of  scatter  were  noted 
in  the  results.  Broutman,  et.  al.  [16]  showed  that  the  experimented, 
data  in  the  tension-compression  quadrant  followed  the  Coulomb-Mohr  theory 
approximately ; 


where  0,  and  a, are  principal  stresses  and  a,,,  and  a,„  are  ultimate 
1  2  *  ut  ue 

tensile  and  coKpressive  strength.  Tiie  modified  maximum  strain  energy 
theory 


(1*3) 


where  V  is  the  poisson's  ratio,  appears  to  serve  as  an  upper  bound. 


CHAPTER  II 


WEAKEST-LINK  STATISTICS  AND  WEIBULL  THEORY 
In  this  Introductory  presentation,  the  concept  of  'weakest-link* 
statistics  for  fracture  [17]  due  to  tensile  stresses  has  been  used  to 
develop  this  theory  of  probability.  Consider  a  volume  V  of  the 
material  which  has  the  property  that  fracture  of  a  sub-volume  v  will 
cause  catastrophic  failure  of  the  whole.  V  is  assumed  to  be  much 
larger  them  v.  For  generality  and  simplicity  elemental  volumes  v^ 
are  assumed  to  be  equal  and  cubic.  Each  of  the  elemental  volumes  will 
fracture  at  certain  stress  O,  characterstic  of  that  element,  and  we 
assume  a  statistical  distribution  of  these  strengths. 

The  probability  density  function  for  the  elemental  volumes  is 
designated  f(a).  It  can  be  defined  in  the  following  way:  from  an 
extremely  large  number  of  elemental  volumes,  the  fraction  will  fail 
between  o  and  0+dc  is  f(a)dc  which  is  the  same  as  the  probability 
that  failure  will  occur  between  O  and  O+dO.  (1 igure  2.1). 


Figure  2.1,  Probability  Density  Function,  f(o) 
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nie  probability  that  failure  will  occur  before  stress 
reached  is: 


g(a) 


a 

j  f(o)do 

O 


a  is 


(2.1) 


The  probability  that  failure  will  occur  before  an  infinite 
stress  is  reached  is  certain,  that  is,  a  probability  of  1.  niis  means 
that  the  function  must  be  normalized  such  that: 


OC 

J  f(o) 

o 


=  g(“) 


-  g(0) 


1 


(2.  2) 


thus  the  probability  that  the  elemental  volumes  will  survive  a 
stress  a  is : 


l'-  g(a) 


(2.3) 


In  the  weakest  link  model.  Figure  2.2  the  volume  V  is  imagined 
composed  of  a  large  number  N  of  elemental  volumes  v,  then  the  number 
of  elements  is  N  =  V/v.  It  is  now  assumed  that  the  failure  of  one 
element  causes  the  failure  of  the  whole.  The  probability  of  failure 
of  the  whole  in  the  stress  range  o  to  o  +  da  due  to  failure  of  one 
element  is  the  probability  of  failure  of  that  element  f(a)da, 

N-1 

multiplied  by  the  probability  that  N-1  elements  survive  {l  -  g(d)}. 

Thus,  probability  of  this  occurrence  is  the  product  of  the 
probabilities : 

N-1 


f(a)da  {  1  -  g(o){ 


(2.4) 


g(a)...  =  Probability  of  Elemental  Failure 
Figure  2.2.  The  WeaJcest-Link  Model 

On  the  other  hand  two  breaks  might  occur  simultaneously  with  the 
probability : 


f(0fda*{l  -  g(o)} 


N-2 


(2.5) 


but  such  simultaneous  events  are  very  rare  and  can  be  neglected.  Thus 
the  probability  density  function  for  the  volume  V  consisting  of  N 


elements  is  proportional  to  ; 


The  elements  that  fall  could  be  any  one  of  the  N  elements ,  thus : 


P(o)  =  Nf(0)  jl  -  g(a)  I  (2,7) 

N 

For  large  N  and  small  g,  the  Poisson  approximation  is  good.  (1  -  g)  = 
-Ng 

e  which  can  be  seen  from  the  Taylor  expansions  of  the  left  and 
right  sides: 


N  N(N-l)g^  N(N-l) (N-2)g^ 

(l-g)  =  1-Ng  +  - - - 

2  3 

2  2  3  3 

-Ng  N  g  N  g 

e  =  1-Ng  + - +  . 

2  3 


(2.8) 


also  for  large  N,  (N-1) =  N  to  make: 

N~1  -Ng 

P(a)  =  Nf (C) (l-g)  =  Nf (0)  e 


(2.9) 


Again  normalization  predicts  certain  failure  for  an  infinite  stress, 
giving: 

OB 

j  P (a) do  =  1  (2.10) 

o 

Now  the  probability  of  failure  function  for  the  voliane  V  coagsrised 
of  N  elements  can  be  written  from  eq.  2.9.  It  has  the  general  shape 


F(o) 


j  P  (a)  do  “  njn 


(a)d0 


o 


o 


(2.11) 


given  by  that  of  Figure  2.3 


Figure  2.3.  Cumulative  Distribution  Function,  Fla) 

Equation  2.9  gives  the  normalized  probability  density  function 
and  equation  2.11  is  the  basic  failure  probability  for  weakest  link 
statistics . 

If  we  consider  the  experimental  strengths  obtained  from  a  series 
of  specimens,  we  will  find  that  there  is  a  definite  relationship 
between  the  probability  that  a  specimen  will  fracture  and  the  stress 
to  which  it  is  subjected.  This  relationship  is  often  called  the 
distribution  f  nction  of  the  strength.  Every  specimen  has  a  unique 
distribution  function.  In  other  words,  the  distribation  function 
is  dependent  upon  the  size  and  design  of  the  spe  i>  i.  Many  different 
functions  have  been  assumed  for  g,  the  failure  probability  for  the 
elements.  The  one  that  is  simple  mathematically,  agrees  with  a  wide 
variety  of  experimental  data,  covering  not  only  material  strengths 
but  many  other  statistical  data  [12J  and  which  appears  physically 
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reasonable  when  the  use  of  it  in  the  above  statistics  is  examined, 
is  the  Weibull  distribucion. 


g(a)  = 


g(a)  =  1 


g(o)  =  0 


j  0<z0^  (2.12) 


Notice  that  when  a  <  O^,  there  is  no  probability  of  failure. 

This  is  the  3-paramoter  distribution.  Some  researchers  assume  =  0 
because  this  leads  to  helpful  simplifications  in  the  application  of 
Weibull  statistics  to  experiments  [is].  A  typical  distribution  curve 
of  strength  [19]  for  =  0  is  shewn  in  the  Figure  2.4.  This  lower 
bound  0^^  on  the  strength  of  the  element  is  also  the  lov/er  bound  for 
the  large  volume  .  must  be  found  for  the  particular  material  and 
for  the  failure  mode  being  considered,  by  experiment.  Because  the 
greatest  value  a  probability  of  failure  function  can  have  for  the 
element  is  can  not  represent  this  function  beyond  a 

certain  value  of  a.  That  value  is  found  from  the  equation: 


g(a  )  = 
max 


(2.13) 


This  will  occur  when  a  =  O  +  o  »  therefore  we  have  put  g(a)  =  1 

max  u  o 

for  a  0^+  Oq  .  Thus  greatest  strength  the 

elemental  volumes  could  have  if  the  Weibull  statistics  were  valid 
over  the  range  0  rZ  0^+  .  However  it  is  only  the  lower  tail 

of  g  that  is  effective  in  *wea>;est-link'  statistics,  and  it  is 


Tensile  Strength  P.S.I. 


Figure  2.4.  A  Weibull  Distribution  of  Tensile  Strengths 

(Source:  Hudson,  J.A.  and  Fairhurst,  Charles, 
"Tensile  Strength,  Weibull’s  Theory  and  A 
General  Statistical  Approach  to  Rock  Failure.") 


reasonable  that  another  simple  power  function  could  be  fitted  to 
the  data  in  this  limited  region.  Using  £g.  2.11  and  2.12  now  the 
probability  of  failure  for  the  volume  V  comprised  of  N  elements 
can  be  written  as : 


FCC) 


(2.14) 


Since  volume  V  is  proportional  to  number  of  elements,  N,  Eq.(17)  can 
be  written  as : 
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A  simple  exan^le  can  show  how  only  the  lower  tail  of  the 


distribution  function,  g,  is  important.  Consider  a  sample  of 

6 

100*100*100  elements.  Thus  N  =  10  .  A  typical  value  for  the  exponent 
m  is  5.  Let  us  calculate  the  probability  that  failure  for  a  volume 
in  the  lower  1/10  of  the  strength  range  of  each  element,  stressed 
unicLxially  and  uniformly.  Inserting  O  =  +  1/10 (Oq)  into  Eq.  2.14 

the  probability  of  failure  function,  using  is; 

,5 

I 

(2.16) 


F(a)  =  1  -  exp  I  -  10“  — 

'10 


=  0.999954 


Thus  the  probability  of  failure  is  almost  1  even  though  the 
threshold  stress  was  exceeded  only  by  0^/10,  Consequently  the  shape 
of  f(0)  and  g{CJ)  beyond  a  +  0^/10  in  of  little  interest. 
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WEIBULL'S  THEORY 

Weibull's  statistical  theory  is  based  on  the  'weakest-link ' 
concept.  The  fact  the  term  'weakest-link'  Is  used  because  the  same 
statistics  apply  to  a  chain  whose  strength  is  described  as  the  strength 
of  its  weakest  link.  In  other  sense,  if  N  seemingly  identical 
specimens  are  broken  in  identically  the  same  way,  it  is  quite  possible 
that  no  two  would  fracture  at  the  same  load.  Instead,  they  can  show 
a  variability  of  100%  or  more;  however,  the  greater  portion  most 
probably  would  fracture  within  a  narrow  range  of  loads,  which  leads 
one  to  believe  that  the  strength  of  a  material  might  lend  itself  to 
statistical  treatment.  To  relate  the  model  of  a  chain  to  that  of  a 
tension  specimen  one  must  imagine  a  multitude  of  fibers,  each  acting 
like  a  chain,  it  can  be  shown  that  the  probability  that  a  model 
specimen  of  this  type  will  fracture  at  a  certain  load  is  a  function 
of  its  volume  [20] . 

As  has  been  pointed  out  above,  a  specimen  may  be  considered  to 
be  made  up  of  meuiy  fibers.  Each  of  these  fibers  can  be  considered  as 
having  its  own  probability  of  withstanding  a  load,  but  the  probability 
for  the  entire  specimen  will  be  the  result  of  considering  the  individual 
probabilities  of  all  the  fibers.  Giving  an  example  of  bending 
specimen,  one-half  of  the  fibers  in  the  test  volume  will  be  in  tension, 
and  the  remaining  fibers  will  be  in  compression.  Weibull's  theory 
cssumes  that  fibers  in  compression  have  a  zero  probability  of  failure. 

As  a  result,  these  compressed  fibers  do  not  contribute  to  the 
probability  for  the  entire  specimen.  Therefore,  an  entire  half  of  the 
specimen  is  neglected  in  Weibull's  consideration  of  the  probability 
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of  a  specimen's  withstanding  a  bending  load.  But  in  considering  the 
probability  of  a  specimen's  withstanding  a  load,  every  fiber  must  be 
considered. 

Using  this  concept  the  Weibull  distribution  function  for 
predicting  the  probability  of  failure  is  expressed  as: 


(2.17) 

V  is  the  volume  of  the  material  under  a  uniaxial  stress  a.  Generally 
0  may  be  a  function  of  the  volume,  in  which  case  F  will  be  a 
function  of  the  external  loading  paramieters . 

The  above  Weibull  distribution  function  describes  the  behavior  of 
many  simple  and  important  load  carrying  elements.  For  a  tension 
member  under  a  constant  stress  state,  the  above  distribution  fxuiction 
can  also  be  expressed  as  follows: 


,  m 

F(o)  =  1  -  exp|-c(o  -  Oy)  }  (2.18) 

1 

where  c  =  -  „ 

^'un(^o^ 

and: 

Vg  =  Gage  volume  of  the  material 
Unit  volume  of  material 
O  =  Applied  uniform  tensile  stress 
F  =  Probability  of  failure 

=  Threshold  stress  below  which  there  is  zero  probability 


of  fracture 
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m  =  Distribution  constant  or  Weibull  modulus 
c  =  Constant 

G^=  A  free  parameter  needed  to  fit  the  function  to  the  data 

0  ,  m,  and  O  (or  c)  are  material  constants, 
u  o 

Parameters  m,  cind  are  material  constants  in  the  true 

sense  and  will  not  change  with  volume  provided  the  Weibull  function 
is  truly  applicable  to  the  material  under  consideration,  eind  if  the 
materials  used  for  small  test  specimens  and  for  the  structural 
components  are  truly  identical.  However,  or  c  is  merely  a 
normalizing  factor  and  is  not  otherwise  related  to  any  physical 
property  of  the  material.  The  value  of  m  is  a  measure  of  the  strength 
varicibility  or  scatter  and  it  helps  in  determining  whether  the  material 
contains  flaws  of  highly  variable  severity  or  similar  geometry  [21]  . 

Estimation  of  Weibull  parameters  from  experimental  data  is  a 
very  important  step  in  the  determination  of  material  characteristics . 

A  large  number  of  methods  are  in  use  to  obtain  Weibull  parameters. 

Heavens  and  Murgatroyd  [22]  in  their  work  used  the  methods  of 
linearization,  direct  curve  fitting,  and  maximum  likehood.  The  problem 
is  complicated  by  the  fact  that  no  method  can  be  proved  to  be  best 
under  all  circumstances.  One  can,  for  example,  estimate  the  mean  of  a 
normal  distribution  by  taking  the  arithmetic  average  or  by  finding  the 
most  frequently  occuring  value  of  the  data.  Davies  [23]  has  suggested 
another  method  for  the  estimation  of  Weibull  parameters.  This  method 
equates  the  mean  variance  and  higher  moments  of  the  distribution  of 
those  of  observed  data.  Another  technique  for  finding  these  parcimeters 
is  based  on  the  moment  generating  method  and  on  rank-order  theory  [24,25], 
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These  results  differ  from  those  obtained  by  Iteration  methods  whereby  the 
sums  of  squares  of  deviations  are  minimized  as  Is  the  case  of  standard 
least-square  techniques.  A  detailed  explanation  of  such  methods  have 
been  given  In  Appendices  A  and  B .  These  methods  can  Involve  large 
truncation  errors  and  can  be  significantly  less  precise  than  conq>arlson 
with  cumulative  totals  (Histogram  Method  for  the  Solution  of  Welbull 
Parameters)  of  experimental  data  (Appendix  C) .  Further,  although  the 
least-square  analysis  method  tends  to  emphasize  deviation  In  the 
experimental  curve  ordinates,  incorporation  of  appropriate  weights  can 
counter  this  effect  while  preserving  the  method's  basic  simplicity. 

Another  problem  with  these  techniques  Is  that  an  assumption  is 
generally  required  for  multiaxial  stress  states  necessitated  by  the 
uniaxial  stress  state  under  which  data  is  taken.  Many  persons  assume 
statistical  independence  in  the  actions  of  principal  direction  stresses, 
i.e.,  the  probability  of  failure  of  an  element  is  equal  to  the  product 
of  probabilities  of  failure  due  to  each  principal  stress  [26,  27].  Under 
three  principal  direction  stresses  the  probability  of  failure  is  given 
by: 

F(o)  =  [F(Oi)]  [Ffo^)]  [FiOj)]  (2.19) 

In  another  theory,  Batdorf  and  Cross  [28]  and  Batdorf  [29,  30] 
assumed  flaws  to  be  penny-shaped  and  neglected  the  shear  on  the  crack 
plane.  Fracture  of  a  crack  was  assximed  to  depend  only  on  the  components 
of  stress  state  normal  to  the  crack  plane.  Freudenthal  [21]  and 
Margetson  [14]  extended  Weibull's  uniaxial  theory  for  stress  assuming 


that  cracks  obey  the  assuaq>tlons  of  Independence  and  shear  insensi¬ 
tivity.  Evans  [31]  developed  a  theory  for  statistical  analysis  of 
fracture  under  multiaxlal  states  of  stresses  based  on  a  critical 


coplaner  strain-energy  release-rate  fracture  criterion. 

Because  a  satisfactory  and  general  treatment  is  not  available  in 
the  literature  particularly  for  the  three-parameter  family,  an  analy¬ 
tical  method  based  on  the  same  least-square  method  using  an  iterative 
programming  scheme  is  developed.  The  objective  is  to  find  those 
values  of  Weibull  parameters  which  make  the  theoretical  curve  fit  the 
experimental  data  bast  for  the  Weibull 's  materials  subjected  to 
various  states  of  stress.  The  methodology  thus  developed  is  illustrated 
by  applying  it  to  two  sets  of  man-made  data  (Chapters  III  and  IV) ,  Addi¬ 
tional  properties  of  Weibull  statistical  distribution  are  given  in 
Appendix  D ,  Including  a  Goodness-of-Fit  analysis  teat  that  supplements 
topical  material  of  this  Chapter. 


CHAPTER  III 


DEVELOPMENT  OP  SIMPLE  TEST  ANALYSIS  CASES 

Weibull's  classical  statistical  theory  includes  two  basic  criteria 
of  failure;  size  and  nonnal  tensile  stress.  Within  the  validity  of 
these  postulates  it  is  capable  of  describing  failure  for  any  type  of 
stress  distribution,  uniform  or  non-uniform,  unicixial  or  polyaxial. 
Failure  in  cui  isotropic  and  homogenous  material  is  more  fully  described 
by  the  three-parameter  family  than  by  the  two- parameter  family.  No 
special  allowance  is  made  for  nonuniformity  of  stress  distribution  other 
than  tiiat  imiilied  in  tiie  iniegration  procedure,  Eq.  3.1  •  Daniel  and 
Weil  [32,  33]  used  standard  specimen  shapes  to  generate  Weibull  statistics 
for  some  typical  cases  of  bending  loading  conditions.  The  three- para- 
met.  '  family  is  developed  here  -  first  for  uniaxial  case. 

Uniaxial  Constant  Stress  Fields  and  the  Determination  of  Weibull 
Parameters 

The  Weibull  distribution  for  a  uniaxial  stress  field  in  a  homo¬ 
genous  isotropic  material  at  a  given  uniform  stress,  0  is  given  bys 


F(a) 


1  - 


exp 


dV 


un 


0 


au<o= 


'V 


u 


OcO, 


(3.1) 
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where:  =  Unit  Volume 

Oy,  Oq,  m  •=  Weibull  paramters  which  are  associated  with  the 
material  and  are  independent  of  siz' 


For  uniaxial  consteint  stress,  o  is  a  principal  stress,  independent  of 
the  volume  so  that  : 


Vg  / 

'  0  -  Ou  ^ 

m 

1  -  exp 

Vun^ 

Oq  y 

(3.2) 


substituting 


c  = 


o 


F(a)  =  1  -  expj-  c(a  -  }Vg  (3.3) 

The  first  solution  procedure  (the  Log-Log  method)  is  based  on 
linear  regression  method  which  converts  the  exponential  type  Weibull 
distribution  function  into  a  straight  line  relationship  by  taking  the  log 
of  the  distribution  function  twice  ,  thus: 

liiln  ( - - \  =  Jln(c)  +mJln(a-a„)  (3.4) 

\l  -  F(a)  / 


This  is  now  a  linear  relationship  allowing  a  linear  least-square 
analysis  to  be  invoked  for  : 


Y 

where  A 
and  B 
If 


A  +  Bx 

J.n(c) 

m 

X  (Oi) 


£n  (Oi  -  Oy) 
inin  I - ^ ^ 

\  1  -  F{0^)  / 


(3.5) 


and  Y^ 


y  (Oi) 


(3.6) 
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i  =  l,  2,  3  n,  where  n  is  the  nuiaber  of  data  points.  Then  A,  B 

are  constants,  which  must  be  determined,  using  the  least-square  crite¬ 
rion.  The  well  known  solution  that  minimizes  the  residual  is; 

(Wi  Xi*)  (rwi  Yi>  -  (ZWi  Xi)  (Ewi  X^y^) 

A  =  - 

D 


(Zwi  Yi  X^)  (Iw^)  -  (Zwi  X^)  (Zwj^  y^) 

D 

(3.7) 


where : 

D  =  (Zw^)  (Zw^  -  (Zw^  X^) ^ 

w.  =  Arbitral-^'  weights  assigned  to  data  points 
All  surr.ations  art  i  =  1,  2,  . .  n. 


In  case  equal  weights  are  desired,  the  quantity  Zw^  is  replaced  by  n. 
The  obvious  difficulty  with  this  analysis  is  that  A  and  B  are  dependent 
on  0Q.  To  determine  the  best  value,  an  iterative  programming 
scheme,  computer  program,  UNIAX.FOR,  was  developed  which  recomputes 
the  residual  sum  of  squares  and  chooses  0^  minimizing  this 
function.  Once  is  known,  m  is  determined  as  the  slope  of  the 
straight  line  expression  and  c  is  found  as  exp (A)  in  Eq.  3.5  . 

A  slight  modification  this  scheme  can  be  used  to  redefine  the  weights 
iteratively  to  emulate  any  curve-fitting  power  law.  This  would 
allow  the  simplicity  of  the  least-square  solution  to  be  invoked 
while  the  weights  are  dynamically  altered  to  (1)  emphasize  the  leading 
distribution  data  or  (2)  de-emphasize  errant  points.  Stamdard 
least-square  techniques  have  the  disadvantage  of  emphasizing  any  point 
away  from  the  curve  by  squaring  that  deviation. 


Comparison  Methods 

As  discussed  earlier,  ocher  techniques  exist  to  determine  Che 
Welbull  parameters  from  the  series  of  data  obtained  from  the 
experimental  tests.  As  an  example,  the  results  from  the  Moment 
Generating  Method  (Appendix  A)  and  the  Newton-Raphson  iteration 
(Least-square)  method  are  compared  with  those  obtained  from  the  Log- 
Log  Method. 

Table  3.1  demonstrates  that  while  a  surprising  discrepancy  can 
exist  in  the  values  of  the  Weibull  parameters  (c,  m,  o^) ,  the  overall 
probabilities  can  appear  very  close  —  at  least  in  the  mid-range.  In 
the  lower  range,  however,  the  probabilities  can  differ  greatly. 
Unfortunately,  this  range  is  of  major  interest  in  dealing  with 
materials  where  the  test  elemental  volume  is  small  in  comparison  to  the 
volume  of  the  structure  being  modeled.  (All  weights  were  assumed  unity 
for  this  example.) 

Pure  Bending  Field  in  a  Rectangular  Section  Bar 

In  the  application  of  brittle  non-metallic  materials,  tests 
utilizing  bending  stress  distributions  to  failure  on  small,  rectangular 
cross-section  bars  are  common  and  generally  less  expensive  than  uniaxial 
tension  tests  —  particularly,  in  the  case  of  three  and  four  point  load¬ 
ing  systems.  In  the  latter  case,  the  central  portion  of  the  bar  is 
subjected  to  a  constant  bending  moment.  Applying  the  Weibull  expression 
[32,  18,  34],  where  the  stress  distribution  is  variable  and  must  be 
integrated  over  the.  volume  in  bending,  the  probability  of  failure  in  the 
pure  bending  section  of  a  uniform  prismatic  beam,  simply-sur  orted  and 


loaded  as  in  Figure  3.1  is: 


125  in 


k .s .i  , 


7.55 

8 . 55 
10.91 
11. 60 
12.44 

13.50 
13.65 
14.00 
14,75 
14.83 
15.96 
16.41 
16.60 
20.00 

21.50 


0.0625 

0.1250 


0 . 2500 
0.3125 
0.3750 
0.4375 
0.5000 
0.5C25 
C.6250 
0 .6875 
0.7500 
0.8125 
0.8750 
0.9  325 


.M.G.Metl'.od 


0.0279 
0.0592 
0.1979 
0.255 
0.5338 


0.4937 

0.5716 

0.5798 

0.6902 

0.7304 

0.7465 

0.9409 

0.9755 


K-R  Method 


0.0362 

0.0653 

0.1913 

0.2453 

0.3208 

0.4284 

0.4444 

C.4S23 

0.5643 

0.5730 

0.6921 

0.7359 

0.7535 

0.9575 

0.9866 


L.L. Method 


0.0618 

0.0970 

0.2258 

0.2757 

0.3429 

0.4353 

0.4489 

0,4808 

0.5497 

0.5571 

0.6583 

0.6966 

0.7122 

0.9190 

0.9631 


WEIBULI,  PARA-METEFi; 

c  (in.  ^  k  ,s  .  i  .  *") 

0.0007079 

0.00002329 

0.00003219 

m  (dimensionless) 

3.0 

4.03 

3.76 

Oy  (k.s.i.) 

4.13 

1.35 

0.026 

o  (k  -s ,1 , ) 

o 

5.61 

14.11 

15.67 

RESIDUES  (dimensionless) 

0.028309 

0.026618 

0.038729 

Table  3.1.  A  Comparison  of  Methods  of  Uci.erniining  Weibull  Parameters 
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Beam  Lioading 


2 


Figure  3.1.  Prismatic  Beam  under  4-Point  Loading  and  Distribution 

of  Extreme  Fiber  Stress 
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Both  support  and  loading  points  are  syiranetrically  positioned  about 
the  beam  center.  The  risk,  of  rupture  in  the  outer  portion  of  the 
specimen  has  been  neglected  in  Eq.(3.8)  but  it  is  easily  Included. 
Fig.  3.1  depicts  this  case  with  its  gage  volume  equal  to  the  pure 
bending  segment.  The  distribution  of  tensile  stresses  in  central 
portion  of  the  test  volume  subjected  to  uniform  bending  is: 

a  2v  _ 2o>57 

=  r  '  (3.9) 

for  k-:'xr^(L-k) 

where  the  extreme  fiber  stress  at  bottom,  for  a  4-point  loading 

is  given  as : 


;  P  =  Fracture  load  (3.10) 

bh* 

This  distribution  follows  the  same  form  of  the  case  of  a  constant 
uniaxial  stress  field  and  is  solved  by  the  same  technique  as  with  the 
uniaxial  constant  stress  field  except  that  different  constants 
result  when  the  log-log  is  taken  in  Eq.  3.8  instead  of  Eq.  312.  The 
programming  scheme  is  given  in  BENDIN.FOR,  Appendix  D,  thus: 


inx-n 


1 

1  -  F(aj^) 


£n(v/2)  -  £n(nH-l)  +  (m+l)  Jln(o^  -  0^) 


-  -  m  h\[a^) 


(3.11) 


In  this  case; 

A  -  Jln(V/2)  -  -  m  in(a^) 

B  =  (m+l)  (3.12) 

with  Oy  and  m  estimated,  is  calculated  from  the  Eq.  (3.12) ■ 
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Pure  Bending  for  a  Circular  Cross-Section  Rod  or  Tube 

This  case  is  similar  to  the  one  discussed  earlier,  except  that 
the  shap>,  of  the  croas-aection  is  circular  .  Davies  [35]  in  his 
work  observed  that  any  minor  defects  incurred  during  manufacture 
along  the  sharp  edges  of  beam  specimens  of  rectangular  cross-section 
that  could  influence  significantly  the  probability  of  failure.  Beam 
specimens  of  circular  cross-section  have  no  sharp  edges  and  are 
therefore  an  attractive  alternative  geometry’. 

Unfortunately  the  use  of  a  simpler  geometry  specimen  leads  to  the 
complication  that  no  closed-form  solution  exists  for  the  probability 
formulation  in  bending. 

The  stress  distrib'Ution  for  tine  case  of  4-point  loading  of  a  rod 
is  : 


a  = 


M  y 


I 


2Pky 

M  trio 


or  for  a  tube : 


o  = 


M  y 
I 


2Fky 

-  ’  k<rx;S:(L-k) 

t  (Ko'"  -Ri*  ) 


(3.13) 


where  P  is  thie  fracture  load. 

For  the  failure  probability  expression  the  integration  is 
conducted  over  the  gage  volume,  v^,  which  is  in  pure  bending.  Starting 
from  the  integral  form  for  tine  risk  of  rupture.  Bn,  Figure  3.2 


Bn 


(3.14) 


Bending  Moment 


Figure  3.2.  4-Point  Loading  for  a  Circular  Cross-Section 

Tube  in  Pure  Bending 


The  limits  of  integration  for  x  in  Eq.  3.14  are  0  and  {L-2k) ;  thus 


Bn 


=  (]>2k)c  f  f  a\ 


dz  dy 


(3.15) 
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1 

where  c  =  - 

The  probability  of  failure  is: 

F(M)  =  1  -  exp (-Bn) 

The  liiaits  of  integration  for  z  and  y  in  Eq.  (3.15)  will  be  different  In 
different  cases.  These  cases  are  considered  as  below. 


(1)  Rod  {Ri=0) 

(a)  0:s:y^<  Ro 

is  the, positive  value  of  y  for  which 
0=0^,  thus: 

Ou  I 

^th=  — 

Parts  of  the  cross-section  for  which  y<y^  inake  no  contribution 
to  Bn.  For  a  fixed  y  such  that  the  limits  of  integration  for 

2  are  obviously; 

-  / (Ro^-  y^)  and  +  /(Ro^-  y^) 


Tfxerefore, 


-  y') 

My  X"" 

—  -  a^j  dz  dy  (3.16) 

-/(Ro^-  y^) 

integration  of  Eq.  (3.16)  with  respect  to  z  gives: 
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Ro 


dy 


(3.17) 


(b)  »  Iri  this  case,  obviously 

Bn  =  0 

(2)  Tube  (Rj:^  0} ;  Three  cases  occur; 

(a)  Ri>y^^>0  ;  The  expression  for  Bn  becorses: 


Ri 


Ro 


//(V  -  "») 


M  y 


Bn  =  (l/-2k)c(  I  J  ^  — T"  -  j  dz  dy  + 

(3.18) 

For  the  first  integral,  the  liadts  of  integration  for  z,  for  fixed  y, 
will  be  from 


/  i  ~  I  z  2 

-  V  (Ro  -  y  )  to  -  y  (Ri  -  y  ) 


and  then  from 


f  2  /  2  2~ 

+  y(Ri-y),  to  +v(Ro-y) 


For  the  second  integral,  the  limits  of  integration  for  z,  for  fixed  y. 


will  be  from 
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Performing  the  integration  gives; 


(b)  Ri:^yj^-cRD 

Following  the  procedure  used  above  gives : 


This  expression  is  identical  to  that  in  the  case  of  a  rod,  Eq.(3.17)> 
In  fact  the  case  of  a  rod  is  obtained  from  this  one  by  setting  Ri  =  0. 
(c)  ;  Here  obviously 


Bn  =  0 


(3.21) 


CHAPTER  IV 


DESIGN  OF  THE  BENDING  EXPERIMENT 


A  major  objective  of  this  contract  was  to  design  tests  and  test 
specimens  that  provided  a  large  amount  of  Information  on  fractlve 
strength  of  brittle  materials  at  low  cost.  As  such,  tests 
such  as  the  uniaxial  tension  test  were  abandoned  due  to  the  need  of 
specially  fabricated  specimens  and  equipment  where  small  strain  to 

failure  ratios  In  ceramics  are  an  Important  design  consideration 
t36J. 


A  four-point  test  fixture  and  tube  specimen  geometry  were 
selected  because: 

1)  A  simple  circular  specimen  geometry  (solid  or  tubular) 
could  be  used  without  excess  concern  of  stress 
concentrations  of  polygonal  section  specimens. 

2)  The  geometry  was  an  inexpensive  one  to  emulate  with 
virtually  all  potential  specimen  materials.  In  many  cases, 
tubes  are  'off-the-shelf’  items. 

3)  Either  smooth  or  rough  surfaces  could  be  tested. 

4)  Without  great  difficulty,  necked  down  specimens  if  needed 
could  be  designed  to  fit  the  same  test  equipment.  The 
need  for  such  designs  would  be  apparent  if  the  frequency 
of  failures  at  the  load  points  is  high. 

5)  Four-point  loading  was  selected  over  three-point  loading 
because  it  provided  a  large  volume  in  the  test  region 
subjected  to  a  known  stress  field. 

Figure  4.1  demonstrates  the  geometry  used  for  the  four-point  load 
test.  The  fixture  components  are: 

A.  Base 

B.  Specimen  resting  on  ground  rollers 

C.  Copper  pads  to  reduce  stress  concentrations  at  bearing  points 
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D*  1.25"-  diameter  alignment  rods  for  stabllzlng  upper  load  frame 

E.  Upper  pivot  arm  with  load  roller  and  specimen  grooved 
rollers 

F.  top  guide  block  above  which  load  Is  applied 

Figure  A. 2  Is  an  assembled  view  with  specimen  In  place.  Figures 
4.3  and  4.4  depict  the  dimensions  and  material  details  of  the  loading 
frame  system  as  finally  employed.  A  load  Is  applied  at  the  top  center 
of  the  top  guide  block  through  a  Tinlus  Olsen  universal  testing  machine. 
This  machine  is  accurately  calibrated  for  load  although  the  rate  of  load 
can  only  be  approximated.  Within  the  machine,  the  fixture  of  Figure  4.2 
is  placed;  loading  rates  were  approximately  40  pounds  per  second. 

During  the  design  and  building  of  the  testing  apparatus,  there  was 
a  question  as  to  what  type  of  alignment  mechanism  should  be  used  to 
locate  the  outer  (lower)  and  the  inner  (upper)  specimen  support  points  so 
that  no  sliding  frictional  constraint  would  be  imposed  during  deflection. 
Stanley  [18]  did  account  for  such  an  effect  as  it  applied  to  a  two-para¬ 
meter  Weibull  test  when  the  coefficient  of  static  friction  between  the 
specimen  material  and  fixture  were  known.  In  this  apparatus,  the  use  of 
rollers,  made  friction  of  minimal  concern.  Ideally,  the  use  of  two 
'formed'  pads  on  rollers  (one  on  the  top  and  one  on  the  bottom)  would 
eliminate  all  non-self  equilibrating  frictional  forces  —  Figure  4.5. 
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Of  further  concern  was  to  determine  a  specimen  contact  design  geometry 
so  that  no  binding  or  specimen  preload  would  be  Inadvertently  Imposed. 

The  problem  was  to  determine  whether  the  apparent  contraction  due  to 
bending  at  the  lower  points  of  support  would  offset  the  apparent  elonga¬ 
tion  due  to  changes  In  slope  at  those  same  points.  An  approximation  to 
this  problem  would  be  to  determine  whether  A  moves  closer  to  or  away  from 
B.  This  would  determine  whether  the  geometry  of  Figure  4.6a  or  4.6b 
would  be  appropriate  for  the  base  of  a  testing  fixture  for  the  specimen 
of  Figure  4.3. 


(a)  (b) 

Figure  4.6.  Possible  Test  Support  Geometries 

To  determine  whether  distance  AB  lengthens  or  shortens,  first  the 
contraction  of  the  neutral  axis  must  be  found;  it  Is  given  by: 


L  L 

AL  =>  r  (ds  -  dx)' 


’J (ds  -  dx)  =  J"  (A+(y' 


)  -  l)dx. 


for 


small  slopes  /y'/«l,  and  A+(y')  ^  ~  1  +  *s(y') 


(4.1) 


(4.2) 


AL 


~  ^ 


o 


(4.3) 


Thus,  for  a  slender  beam  a  contraction  of  any  given  span  due  to  bending 

However,  the  change  in  slope  of  a  beam  imder  load 


would  always  occur. 


1 
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tends  to  place  the  lower  points  of  the  beam  in  like  curvature  as 
indicated  below  in  Figure  4.7. 


y  (a)  y  (b) 

Figure  4.7.  Bending  Causing  Outward  Movement  of  Beam  Lower  Surface 

Thus,  the  expansion  or  elongation  between  a  and  b  is: 


Li 

ymax  “  yUa)]'  hj 


(y 


*)  ^dx 


(4.4) 


at  the  bottom  of  the  beam  [ymaxl  the  total  offset,  Ro,  wherever 

contact  occurs. 

In  the  case  of  a  hollow  rod  in  pure  bending  for  its  mid-length  L,  we 


have: 


and 


/ 


\  I  (y')  dx 


24e"i^ 


ymaxty‘^^>  = 


W 

El 


RqMqL  _  >  Q 

El  24Eh'  ^ 


(4.5a) 

(4.5b) 


and  one  must  check  for' 


(4.6) 
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For  the  center  span  of  a  bar  specimen  In  four-point  loading,  there 
Is  no  need  to  perform  a  calculation  to  determine  the  geometry  of  the  load¬ 
ing  support  as  both  the  effect  of  slope  and  shortening  of  the  neutral 
axis  dictate  that  the  Figure  4.6a  Is  appropriate  for  the  overall  shortening. 

Four-point  loading,  the  above  equation  Is  an  approximation  to  the 
change  In  distance  between  C  and  D  of  Figure  4.8  below  provided  H  ->■  I 
in  length. 


Figure  4.8.  Bending  Geometry  and  Notation 


Thus,  if  il  -►  £,  then  from  4.5  and  4,6: 

LMo  ~  Mo^L*  . 

El  ^  ■  24  gziz 

(4.7) 

Here  y  Is  the  distance  from  neutral  axis  to  point  of  support  and 


Mo  =  Constant 


PCl-I) 

4 


(4.8) 


The  exact  solution  is  obtained  by  substitution  of  the  moment  distri¬ 


bution  along  the  entire  beam  into  (4.4)  and  Integrating.  In  this  case: 
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El  y'(x)  ■  ^  [<x+m>^  -  <x+m>^  -  <x-m>  -  (m-m)  (m  +  m)  ] 


where  m  ••  il/2,  m  •>  %ll  and  <x  -  Xq> 


Ux  -  Xq)  = 

V 


if 

if  x<x„ 


(4.9) 


Substitution  of  EIy"(x)  derived  from  (4.9)  and  integration  of  each  of 
the  resulting  expressions  leads  to: 

AL  -  -  SOm^m^  +80™“  -  32m®]  (4.10) 

for  the  beam  in  4-point  loading. 

The  slopes  at  the  outer  supports  are  given  by: 

y'(-m)  =  -  (m^  -  S*) 

and  (4.11) 


Thus,  if  the  specimen  is  supported  at  an  offset  of  R  from  the  neutral 
axis.  Figure  4.4a  and  b  ,  then  the  restraining  geometry  Figure  4a  should 
be  used  if  X>o,  otherwise  4b  should  be  used  where: 


X  =  (m^  -  m^)Ro  -  TTV  (32m*  -  80  m*m^  +  80mS  -  32m  )  (4.12) 

4EI  oO 


Thus,  the  specimen  loading  design  is  chosen  with  attention  to 
several  important  details.  By  the  use  of  bearings  at  the  load  points,  the 
effect  of  friction  is  minimized  and  proper  definition  of  the  loading  frame 
geometry  (Figure  4.4a  and  b)  eliminates  a  chance  of  superposition  of 


compression  in  the  test.  Finally,  there  is  a  complicated  stress  field  in 
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vicinity  of  the  four  loading  areas.  Machined  surfaces  on  the  test  fixture 
at  such  points  and  soft  copper  pads  reduced  this  problem,  particularly,  near 
the  two  Inner  contact  points  adjoining  the  gage  length. 
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EXPERIMENTAL  DETAILS 


Specimens  and  Material 


All  specimens  used  in  the  experimental  work  were  tubes  of  circular 


cross-section  made  of  99.8%  Aluminum  Oxide,  AJI2O2,  (998  Alumina,  McDannel 
Refractory  Procelain  Company,  Beaver  Falls,  Pennsylvania,  15010)  or  of 
mullite  (MV33)  from  the  same  source.  The  test  schematic  appears  below 
in  Figure  4.9.  The  tube  geometries  for  both  materials  are  similar  with 
nominal  values  listed  below  Figure  4.9,  and  with  eccentricities  and 
standard  deviations  given  in  Table  4.1. 


Testing  Procedures 

The  general  procedure  followed  was  ASTM  Standards  for  rods,  1977, 


Part  17,  pages  104-11. 


P 


Figure  4.9.  Testing  Procedure  for  Bending  Specimen 
Observations : 

Referring  to  Figure  4.9,  general  dimensions  are  in  inches: 
L  =10 

K  =  3.6 


(L-2k) 


2.8 
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Ro 

■ 

.4999 

(Alumina) 

Ri 

S 

.3091 

(Alumina) 

Vg 

- 

.6789 

(Alumina) 

(under  tension) 

Ro 

- 

.5132 

(Mullite) 

Rl 

- 

.3661 

(Mullite) 

Vg 

.5689 

(Mullite) 

(under  tension) 

998  Alumina 

Mullite  (MV33) 

Mean  and  SD  of  Maximum 

OD 

1.0010+  .0021 

1.0290+  .0065 

Mean  and  SD  of  Direction 
+  to  Max.  OD  Direction 

.9985+  .0027 

1.0236+  .0082 

Mean  and  SD  of  Minimum 

ID 

.6156+  .0203 

.7276+  .0087 

Mean  and  SD  of  Direction 
+  to  Min.  ID  Direction 

.6209+  .0234 

.7367+  .0074 

Table  4.1.  Specimen  Means  and  Standard  Deviations  In  Inches 
for  Inner  and  Outer  Diameters 


N  =  50  (Alumina) 


N  =  50  mullite 
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Analysis  of  Test  Results; 

The  technique  used  In  the  computer  program  HOLLOW. FOR  and  DIN. FOR, 

Is  an  Iteration  method  where  the  values  of  Welbull  parameters,  c,  m,  and 
are  found  by  the  minimization  of  residual  errors  In  the  experiment¬ 
ally  determined  fracture  probabilities  compared  to  £q.  4.10  through  4.12, 
the  probability  of  fracture  F(a)  Is  calculated  by  repeated  applications 
of  Simpson's  rule  [37]  to  the  Integrals  until  a  truncation  error  condi¬ 
tion  Is  satisfied.  Residual  error  as  calculated  by  the  least-square 
method  Is  below: 


2 

Res(c,m,CL)  *  w^{F(a)  -  F(a)}  (4.13) 


\diere  F(a)  =  [Mean  Rank  Method  for 

n+1 

This  process  perturbs  values  of  c, 

minimized.  This  gives  the  final  values 

2-Parameter  Family 

c  =  0.403E-27  In. lbs."™ 
,391E-46  N"™ 

ra  =6.20  (dimensionless) 

=  0 . 0  p . s . i . /MPa 

o  =  26,2i2  p.s.i. 

°  180.73  MPa 

Res  =  0.02434 


the  experimental  fracture  loads] . 

m,  and  o\,  until  the  Eq.  4.13  is 

of  Welbull  parameters  for  Alumina. 

3-Paraineter  Family 

c  =  0.403E-27  in.^™"^lbs."™ 
.391E-46  M^“"^  N"™ 

m  =6.20  (dimensionless) 

=  0.0  p.s.i. /MPa 

a  =  26,212  p.s.i. 

°  180.73  MPa 

Res  =  0.02434 


Table  4.2;Weibull  Parameters  for  Alumina  998 
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The  experimental  results  are  given  In  Table  4.3.  This  table  depicts  | 

failure  probability  estimates  via  the  two-  and  three-parameter  Welbull  I 

distributions  In  comparison  to  the  mean  rank  experimental  results.  Table  4.4 

depicts  those  specimens  that  failed  abnormally  and  were  not  used  for  data  | 

i 

I 

collection.  The  residual  Is  minimized  with  the  two  parameter  family,  thus, 

Table  4.2  entries  are  Identical.  Another  Illustration  Is  given  where  the 

I 

Welbull  parameters  of  Table  4.2  are  assumed  for  a  one  cubic  Inch  specimen  { 

under  uniaxial  stress  states.  Thtis,  a  one  cubic  Inch  specimen  exhibits  a 
50%  chance  of  fracture  if  stressed  uniaxially  and  uniformly  to  24.6  ksi, 

i 

i 


I 


Specimen 

Designation 

Fracture 

Load 

P* 

(lbs.  ) 

Moment 

M 

PK/2* 

( lb-in ) 

Rank 

i 

F(M) 

i/n+l 

Fract'ire  Probability  at  j 
Failure  Moment,  M  1 

F(M)  1 

2-Parameter 

3-Parameter 

A2 

1010 

1818 

1 

0.04762 

0.024402 

0.024402 

i  AlO 

1210 

2178 

2 

0.09524 

0.074186 

0.074186 

A22 

1365 

2457 

3 

0.14286 

0.15019 

0.15019 

A23 

1395 

2511 

4 

0.19047 

0.16992 

0.16992 

A5 

1475 

2655 

y 

0.23809 

0.23168 

0.23168 

A4 

1510 

2718 

6 

0.28571 

0.26273 

0.26273  ! 

A9 

1610 

2898 

7 

0.33333 

0.36467 

0.36467 

A17 

1620 

2916 

8 

0.38095 

0.37585 

0.37585 

A21 

1730 

3114 

9 

0.42857 

0.50773 

0.50773  i 

A18 

1745 

3141 

10 

0.47619 

0.52655 

0.52655  1 

All 

1750 

3150 

11 

0.52380 

0.53284 

0.53284  i 

A1 

1785 

3213 

12 

0. 57143 

0.57705 

0.57705  { 

A16 

1865 

3357 

13 

0.61905 

0.67674 

0 . 67674 

A14 

1875 

5375 

14 

0 . 66666 

0.68881 

0.68881 

A25 

1885 

3393 

15 

0.71428 

0.70076 

0.70076 

A6 

1905 

3429 

16 

0.76190 

0.72420 

0.72420 

A12 

1915 

3447 

17 

0.80952 

0.7358 

0.7358 

A3 

2005 

3609 

18 

0.85714 

0.82948 

0.82948 

A15 

2140 

3852 

19 

0.90476 

0.929235 

0.929235 

A19 

2250 

4050 

20 

0.95238 

0.97311 

0.97311 

Table  4.3  Experimental  Results  for  4-Polnt  Bending  Test  for 
A  Circular  Cross-Section  Tube  of  Alumina  998 


*  English  system  used  since  equipment  Is  calibrated  In  pounds. 


Specimen 

Designation 

Fracture 

Load  (lbs. ) 

Why  Rejected 

A7 

1770 

Failed  out  of  gage  length 

A8 

1690 

Testing  Machine  Failed  Then 
Pulsed  Causing  Shock 

A13 

1640 

Failed  Under  Support 

A20 

1440 

Failed  Under  Support 

A24 

1595 

Failed  Under  Support 

Table  4*4.  Experimental  Results  For  Rejected  Alumina  993  Bend  Test  Specimens 

Figure  4.10  demonstrates  that  the  Alumina  fracture  curve  is  a  relatively 
broad  one  typical  of  experiments  in  which  care  is  taken  in  the  experimental 
design.  Frequently  attempts  to  load  brittle  specimens  in  pure  tension 
result  in  an  even  broader  app.crent  curve  than  the  one  of  Figure  4.10  because 
there  are  inherent  small  misalignments  in  the  experimental  equipment.  These 
aberrations  produce  stresses  that  reconfigure  material  exhibiting  plasticity. 
Brittle  materals  are  subjected  to  stresses  which  only  degrade  by  fracture 
at  a  higher  level  than  moan  or  macroscopic  value  sensed  by  the  recording 
equipment.  The  Weibull  modulus  is  a  recognized  measure  of  the  combined 
uniformity  of  the  experimental  design  and  specimen  strength  predictability. 
For  this  material  and  test,  a  value  of  6.2  is  higher  than  typical. 


aamm  jo  ixniavaosj 
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The  "^sults  of  the  mulllte  bend  tests  are  given  in  Table  4.5  with 
Table  4.6  representing  the  22  specimens  used  to  determine  the  para¬ 
meters  of  Table  4.5.  Table  4.7  Is  a  list  of  rejected  specimens  of 
which  only  two  occurred. 

The  mean  load  at  failure  of  the  mulllte  at  843  poimds  Is  signifi¬ 
cantly  below  the  corresponding  load  of  alumina  at  1,702  pounds.  It 
should  be  noted  that  even  though  the  mulllte  Is  about  one-half  the 
strength  of  alumina,  the  dispersion  to  strength  ratio  of  both  materials 
Is  about  the  same. 

As  with  alumina,  the  mulllte  tests  were  used  to  generate  Uelbull 
parameter  sets  to  produce  the  values  below: 


2-Parameter  Family 

3-Parameter  Family 

c 

2m- 3  -m 

.725E-20  In.  lbs. 

c 

. 

.725E-20  ln.^"'-*lbs."“ 

,  2m-*  -m 

(.6874E-34M  N  ) 

—in 

(.6874E-34M  N  ) 

m 

= 

4.90 

m 

ss 

4.90 

Ou 

= 

0.0  p.s.l./M  Pa 

Ou 

B 

0.0  p.s.l./M  Pa 

= 

12,887  p.s.l. 

(88.856  M  Pa  ) 

Oq 

B 

12,887  p.s.l. 

(  88.856  M  Pa  ) 

Ses 

s 

.0544 

Res 

s 

.0544 

Table  4.5:  Welbull  Parameters  for  Mulllte  (MV33) 

Once  again  the  two  parameter  Welbull  distribution  was  found  to  minimize 
the  results  as  well  as  any  three-parameter  distribution. 


Specimen 

Designation 

Fracture 

Load 

P* 

(lbs. ) 

Moment 

PK/2* 

( lb-in) 

Rank 

i 

F(M) 

i/n>l 

Fracture  Probability 
Failure  Moment,  M 

F(M) 

2-Parameter 

3-Parameter 

M23 

530 

954 

1 

0.0526 

0.0742 

0.0742 

M24 

550 

990 

2 

0.1053 

0.0882 

0.0742 

M5 

595 

1071 

3 

0.1579 

0.1270 

0.0742 

M19 

710 

1278 

4 

0.2105 

0.2764 

0.0742 

M21 

740 

1332 

5 

0.2632 

0.3272 

0.0742 

Ml 

770 

1386 

6 

0.3158 

0.3824 

0.0742 

MIO 

800 

1440 

7 

0.3684 

0.4407 

0.0742 

M9 

810 

1458 

3 

0.4211 

0.4608 

0.0742 

M14 

825 

1485 

9 

0.4737 

0.4912 

0.0742 

M17 

830 

1494 

10 

* 

0.5263 

0.5014 

0.0742 

M22 

830 

1494 

10 

* 

* 

0.5263 

0.5014 

0.0742 

M18 

850 

1530 

11 

0.5739 

0.5426 

0.0742 

M2 

860 

1548 

12 

* 

0.6316 

0.5632 

0.0742 

M4 

860 

1548 

12 

* 

* 

0.6316 

0.5632 

0.0742 

MS 

860 

1548 

12 

* 

* 

0.6316 

0.5632 

0.0742 

M13 

875 

1575 

13 

0.6842 

0. 5940 

0.0742 

M16 

970 

1746 

14 

0.7368 

0.7756 

0.0742 

Mil 

1000 

1800 

15 

* 

0.7894 

0.8236 

0.0742 

M12 

1000 

1800 

15 

* 

X 

0.7894 

0.8236 

0.0742 

M7 

1050 

1890 

16 

0.8421 

0.8896 

0.0742 

M6 

1100 

1980 

17 

0.8947 

0.9372 

0.0742 

M3 

1130 

2034 

18 

0.9474 

0.9575 

0.9575 

*  English  system  used  since  equipment  is  calibrated  in  pounds. 
**  Multiple  entries  treated  with  proportionately  higher  weights. 


Table  4.6:  Experimental  Results  for  4-Point  Bending  Test  for  a 
Circular  Cross-Section  Tube  of  Mullite  MV33 


I 

I 


Specimen 

Maximum 

Designation 

Load 

Remarks 

M14 

840 

Fractured  out  of  gage  length 

M20 

1050 

Fractured  out  of  gage  length 

Table  4.7:  Experimental  Results  for  Rejected  Mullite  MV33 
Bend  Test  Specimens 


The  graph  for  the  Mullite  parameters  corresponding  to  Figure  4.10 
for  Alumina  appear  as  Figure  4.11.  The  slgnlflcatly  lower  strength 
of  mullite  can  easily  be  demonstrated  between  the  two  materials  when 
Figures  4.10  and  4.11  are  compared. 


a^mivi  io  Ainiavso»d 
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Examination  of  Fractured  Specimens 

Some  alumina  specimens.  Figure  4.12,  failed  through  a  single  crack 
wholly  within  the  gage  length  while  no  mulli'-e  tubes  failed  wholly  in  a 
single  crack  pattern. 

The  least  complex  crack  surface  of  a  mullite  specimen  was  compound 
where  1/2  inch  wide  sector  was  separated.  Figure  4.13,  even  then, 
failure  was  wholly  within  the  gage  length. 

The  mullite  failed  in  every  case  with  more  splintering  and  cnip- 
ping  than  the  alumina.  Thus,  all  mullite  specimens  were  taped  for 
safety  and  retrieval  of  splinters.  IVhile  Mullite  fractures  were  more 
complex,  the  alumina  specimens  exhibited  the  most  unpredictable 
failure  geometries;  they  ranged  from  the  simpie  fracture  of  Figure  4.12 
to  fractures  where  a  variety  of  compound  shapes  occurred  designated  by 
multiple  tension  compount  fracture  (Figure  4.14),  multiple  compression 
compound  fracture  (Figure  4.15),  and  two  section  fracture  with  a 
compression  axial  split  (Figure  4.16).  In  no  cases,  however,  did 
noticeable  spalling  or  chipping  of  other  than  material  powder  occur. 

On  the  other  hand,  the  mullite  samples  all  exhibited  a  smaller 
variety  of  fracture  patterns  vith  apparent  single  fracture  initiation 
locations  on  the  tension  side  with  eminating  ray-like  cracks  to^i’ard 
the  compression  side  similar  to  Figure  4.15  in  the  alumina  tests. 

These  mullite  figures  are  listed  in  increasing  order  of  complexity, 
Figures  4.17,  4.15,  and  4.19. 

In  Figures  4.13  and  4.18,  there  is  clear  evidence  of  mullite  spall¬ 
ing  in  the  compression  regions.  Of  course,  strain  energy  densities 
were  highest  on  the  upper  and  lower  regions  of  the  specimens. 
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While  it  cannot  be  seen  clearly,  many  of  the  alumina  specimens 
have  crack  bifurcations  in  regions  where  there  is  an  abrupt  change  in 
crack  path.  For  example,  in  Figures  4.14  and  4.16  where  the  'major' 
crack  path  has  propagated  extensively,  the  crack  has  not  opened  through¬ 
out  the  material.  In  strong  light  specimens  4.14  and  4.16  reveal  these 
bifurication  points  and  hidden  cracks  near  the  neutral  axis  plane  again 
where  the  visible  cracks  abruptly  change  direction.  No  such  phenomenon 
was  detected  in  any  of  the  mullite  specimens  althougVi  they  were  more 
opaque . 


Figure  4.12.  Alumina  Speclracn  bl  After  Fracture 


Flrure  4.1J:  MulJite  Specimen  M2i  After  Fracture 


A17 


Fi  i'u: 

re  H . lb 

:  liv'o-Sec 

tion  Fracture  with  a 

Comp  1 

resslo 

Split  i 

n  Alurcina  Specimen  A 

1  7 

Mullitfc  Bend  Specimen  M19  Multiple  Fracture 


KJ  V 


Figure  4.19:  Mullite  Bend  Specimen  M20  Kultiple  Fracture 


CONCLUSION 
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Vfhlle  It  was  originally  expected  that  a  three-parameter  Welbull 
family  would  be  required  to  describe  adequately  the  probability  of 
fracture  in  both  Alumina  and  Mulllte,  it  was  found  that  the  2-parameter 
family  produced  the  same  minimum  residual  sdm  of  squares.  As  a  result, 
the  parameter  family  was  disregarded  in  favor  of  the  simpler  distribu¬ 
tion  for  all  bending  tests. 

The  mean  failure  load  of  the  Alumina  tubes  was  1702  lbs  in  four-point 
loading  vs  843  lbs  for  Mullite.  That  is  equivalent  to  a  stress  level  of 
170  MPa  (24.6  ksl)  for  Alumina  and  83.4  MPa  (12.1  ksi)  for  Mullite  in 
a  one  (1)  cu..  in.  specimen  of  uniform  stress  at  a  50%  chance  of  failure. 

While  it  is  appropriate  to  employ  specimens  designed  to  be  operated 
on  precision  aligned  axial  tensil  testing  equipment,  the  use  of  these 
machines  and  costly  specimens  have  no  advantages  over  properly  designed 
four-point  loading  tests.  Few  specimens  demonstrated  cracks  or  crushing 
near  the  points  of  supports. 

While  the  mathematics  required  in  the  minimization  of  residuals  for  a 
least  square  method  solution  is  not  closed-form  in  the  case  of  rods  or 
tubes,  the  tt.^.ting  simplicity  compensates  for  the  minor  problem.  It  is, 
however,  true  that  the  computer  model  producing  the  Weibull  paramenters 
could  be  made  to  operate  more  efficiently  and  automatically. 

In  order  that  various  hypotheses  can  be  tested  in  dealing  with 
failure  from  multiaxial  stress  states,  several  stress  tensor  ratios 
must  be  employed.  The  subsequent  or  follower  contract  report 
N00019-79-PR-RL218,  will  address  parallel  work  in  torsion. 
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APPENDIX  A 

RELIABILITY  PREDICTIONS  BASED  ON  TEST  SPECIMENS 

In  the  circumstances,  when  it  is  not  practical  to  investigate  the 
integrity  of  a  component  through  tests  on  small  models ,  the  reliability 
predictions  can  be  based  on  data  obtained  from  test  specimens.  This 
method  requires  the  development  of  an  approach,  where  a  combined-fracture 
stress  theory  is  needed  to  relate  the  behavior  of  a  unit  volume  under 
a  complicated  stress  state  to  the  behavior  of  the  test  specimen  under  a 
controlled  stress  state.  Development  of  such  an  approach  Is  far  from 
complete;  however,  based  on  a  series  model,  Barnett  et.al.  (26)  proposed 
a  simple  theory  which  is  often  used  to  predict  behavior  of  a  volume 
subject  to  a  general  stress  field. 

If  a^,  a^,  and  Oj  are  principal  stresses  acting  on  a  uniformly 
stressed  unit  volume  ,  the  reliability  of  the  volume  is: 

1  -  F(a)  =  [l  -  F{Oj)]  [l  -  F(a2)]  [l  -  FCOj)]  (A.l) 


where ; 

1  -  F(a)  =  reliability  of  the  volume 

F(a)  =  The  fracture  probability  of  a  unit  volume  under  a  pure 
tensile  stress  o 

Eq.  (1.1)  can  be  used  for  any  uniformly  stressed  basic  unit  for 
which  F(a)  has  been  established,  including,  for  example,  infinitesimal 
volumes.  It  is  however  an  assumption  of  independence  of  effect  due  to 
the  stress  tensor  components  actions. 

The  simplest  way  to  determine  the  reliability-strength  trade-off 
is  to  test  full-scale  prototypes,  where  neither  a  stress  nor  strength 
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analysis  is  required,  but  the  drawback  is  expense.  The  testing  of 
scale  models  costs  less  and  still  avoids  the  need  for  stress  analysis. 
However  this  technique  requires  a  knowledge  of  strength-volume  relation¬ 
ship.  Since  neither  of  these  approaches  provides  much  basic  information 
about  the  material  being  tested,  therefore  subsequent  designs  with  the 
same  material  must  be  developed  from  scratch. 

The  most  sophisticated  approach  -  the  statistical  theory  of  which 

is  explained  above  -  requires  not  only  data  from  test  specimen 
\ 

requires  a  detailed  stress  analysis.  Small,  relatively  inexpensive 
specimens  are  tested  to  determine  the  reliability  characteristics  of 
small  volumes  which  further  leads  to  predict  prototype  behavior.  The 
entire  analysis  requires  six  steps. 

1.  Obtain  the  tensile-strength  distribution  F  (a)  using  a 
tension  specimen  with  gage  volume  Vg . 

2.  Perform  a  complete  stress  analysis  of  the  component. 

3.  Divide  the  component  into  n  convenient  volumes,  . ,V^. 

Each  volume  should  be  sufficiently  small  so  that  it  contains 
approximately  a  constant  stress  state. 


4.  Determine  the  'worst*  stress  condition  in  each  volume  and 
assume  that  the  corresponding  principal  stresses  Oj,  Ojr  and  Oj 


act  uniformly  throughout  the  volume. 


5. 


Determine  the  reliability  of  each  volume  (l-F^)^  by  first 
finding  the  reliability  of  the  gage  volume  under  the  principal 


stresses  through  the  application  of  Eq, 


(1.1)  ; 


(A.  2) 


Then  scale  gage-volume  reliability  to  find  d-F^)  by  using 


the  Eq.  (1.3)  . 


v./v 

J  9 


(A. 3) 


6.  Use  Eq.  (1.4)  which  is  based  on  weakest-link  hypothesis,  to 
establish  the  reliability  of  the  entire  structure,  (l-F(a) 
from  the  reliability  of  the  volumes  . 


l-F(a)  =  (1-F.)(1-F-) . “  (1-Fd 

3=1 


(A. 4) 


l-F(a)  =  I  (1~F,  )  . 

3=1  V  3 


(A. 5) 


*  Eq.  (1.3)  is  based  on  the  extreme  value  statistics,  which  furnish 
with  an  important  necessary  condition  for  a  series  material  ttat 
does  not  require  the  specific  form,  of  the  distribution  f\inction  nor 
a  knowledge  of  the  combined-stress  theory  appropriate  for  the 
material.  Specifically,  when  the  loading  and  geometry  of  two  different 
size  components  are  similar,  their  distribution  function  F(0)  must 
scale  according  to  the  following  relationship  when  the  material  obeys 
the  series  model . 

,  v,/v, 

1  -  =  1  -  F,(Oj) 

where  is  the  fracture  probability  of  the  i^  structure,  is  the 
volume  of  the  i*^  structure. 
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APPENDIX  B 

Newton-Raphson  Iteration  Method  for  Matrix  Equations  of 
Welbull  Parameters 


The  residue  equation  is: 


2 

R(c,m,Oy)  -  -  I  1  -  «xp{-c(a  -  Oy)  }|]  (B.l) 


R  =  Residue  as  function  c,  a,  and  a 
where  “ 

F^*»  Data  Points 


Simplification  of  Eq.  (2.1)  takes  the  form  of: 


-  .  ~ui  ui  =  , 

R  -  -2  Z  I  F^e  +  e  +  F^* 


(3.2) 


where  F^  =  (1  -  F^) 
Ui  =  cCa^  - 


Minimization  of  residue  R,  requires  that  partial  derivatives 
uf  Eq.  (2.9)  with  respect  to  c,  m,  and  0^  be  zero. 


—  *  2  Z  Wi  I  F,e  -  e 


'Ui 


(Oi  -  o^) 


(B.3) 


6r  ^  I  =  -Ui  -zui  )  ® 

^  =  2  .Zi  Wi  j  F^e  -  e  [  c  in(Oi  -  0^) (a^  -  O^) 


—  »  -2  Z  w.  {  F.e  -  e 

6au  i=l  ^  ^ 


c  m  (a^  -  0^) 


(B.5) 
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Again,  a  Newton-Raphaon  solution  of  c,  m,  and  a^  requires  that 
the  partial  derivatives  of  Eqs.  ilB.3),  and  (B.5)  with  respect 

to  c,  m,  and  should  be  satisfied.  The  following  Is  the  desired 
matrix  expressed  In  a  symbolic  form. 


R,C 

R,cc  R,cm 

\ 

6c 

\ 

\ 

\ 

\ 

R,m 

=  - 

R.mm  R.mD 

6m 

\ 

\ 

o  \ 

\ 

Or 

Q 

c 

\ 

\ 

L.  -J 
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where  in  condensed  form 


-  c(a.  -  a^)“ 


R,„.  =  2  Z  w,  F,e  -  e 


»  -Ui  -2Uj^  I  <Sui 


j-  (  i 


1  .—i 

’  6v. 


=  -Ui  -2Ui  j 

^'Vjvj,  =  2  Z  Wi  j  FiO  -  e  [ 


6v.6v 
J  k 


^  -2Ui  «=  -Ui  >  6ui  6ui 

E  “ij2«  -^1. 


(B.7) 


Vj  =  one  of  {c,  m,  0^}  ,  j  =  1,  2,  3. 

Solution  of  this  matrix  gives  the  corrections  c,  m,  emd  to 
initial  guesses  of  c,  m,  and  Refined  values  of  W'eibull  parameters 

are  obtained  by  adding  these  correction  values  to  the  initial  guesses 
and  repeating  until  the  values  converge  (computer  program  NEWTON. FOR) . 

Unfortunately,  this  method  is  highly  unstable  and  frequently 
only  converges  if  the  initial  values  are  within  ±  5%  of  the  correct 
values.  The  problem,  therefore,  is  to  produce  very  accurate  initial 
guesses  easily  or  abandon  this  technique  as  a  less  useful  one.  One 
mechcinism  albeit  not  entirely  successful  of  selecting  initial  values 
is  listed  below; 

Initial  Prediction  ^uld  Stepping  Method 

In  this  method  a  set  of  three  points  o^,  CJg,  and  is  chosen 
for  the  analysis  of  Weibull  parameters.  We  let 
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where  n  is  a  multiple  of  4. 


Figure  B.l.  Selection  of  Three  Data  Points 


This  leads  to  a  technically  simple  approximation  for  c,  m. 
Expressing  the  Weihull  distribution  for  point  ,  we  obtain 

=  1  -  ejq)  |-c(a^  -  a^)  | 


and  similar 
Taking 


expressions  for  points  Og, 
In  In  we  have: 


and 


£n  c  +  m  -  a^)  “ 

similcirly  for  points  a„,  and  a_  we  have: 

£.n  c  +  m  ^nCOg  -  a^)  »  Kg 
Jin  c  +  m  Jin  (a  -  Oy)  =  K 


(B.8) 


and  Oy. 

(B.9) 

(B.IO) 

(B.ll) 

(B.12) 
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where 


=  £niln 


1  -  F, 


and  similarly  for  Kg  and  K^. 


Simplification  of  Eqs.  (2.10),  (2.11),  and  (2.12)  gives: 


by  eliminating  m  we  obtain: 


(b.13) 


(B.14) 


(B.15) 


where 

Kb  -  Kc 

y  =  - 

Ka  -  Kb 


Thus  the  problem  is  reduced  to  finding  Oy  of  Eq.  (B.15).  The 
Newton-Raphson  iteration  method  for  the  solution  of  parameter  requires; 


F(au) 

'^“n  “  °“o  '  F'(aj 


(B.16) 


where 


X 


New  value  of  after  each  iteration 
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■  Preceding  value  of  C7y  used  for  reiteration 
F'  (a„)  ■  First  derivative  of  F(o  ) 


Expansion  of  Eg.  (2.16) takes  the  form  of: 


Y  In 


°Uo  - 


+  iln 


V  -  ‘=’b 


r 

V 

\-‘^B 

“  ^B 

■ 

V'b 

T 

. 

Q 

C 

O 

1 

T 

%  -  <^C 

(B.177 


solution  of  the  above  expression  gives  the  value  of  by  iteration. 

Once  Oq  is  determined,  the  values  of  c  and  m  can  be  found  from  Eqs. 

(B.IO)  and  (B.135. 

Unfortunately,  it  can  be  shown  that  solution  of  Eg.  (B,17)  can 
lead  to  Oy  which  generates  complex  c  and  m  even  for  O  «  0^,  Og,  or  Og,  if 
real  probabilities  are  generated.  When  real  c  and  m  are  found,  it  is 
still  often  that  the  Newton-Raphson  process  does  not  converge,  further 
Improvement  in  solution  technique  is  needed.  The  iterative  technique 
is  included  as  a  computer  program  in  Appendix  E,  program  5,  NEWTON. FOR. 
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a  *  Preceding  value  of  0„  used  for  reiteration 
“o 

F'COy)  =  First  derivative  of  F(o^) 

Expansion  of  Eq.  (2.16) takes  the  form  of: 


O 


“n 


j  ■  °A 

1 

%  -  ' 

1  1 

%  “ 

1 

V 

°A-°B 

1 

a  -  a 

B 

T 

- 

■  °A 

T 

J 

(B.17'J 


solution  of  the  above  expression  gives  the  value  of  Oy  by  iteration. 

Once  Oy  is  determined,  the  values  of  c  and  m  can  be  found  from  Egs. 

(B.IO)  and  {B.13  . 

Unfortunately,  it  can  be  shovm  that  solution  of  Eq.  (B.J7)  can 
lead  to  Cy  which  generates  complex  c  and  m  even  for  0  =  0^,  Og,  or  Oq,  if 
real  probabilities  are  generated.  When  real  c  and  zn  are  found,  it  is 
still  often  that  the  Newton-Raphson  process  does  not  converge,  further 
inqzrovement  in  solution  technique  is  needed.  The  Iterative  technique 
is  included  as  a  computer  program  in  Appendix  E,  program  5,  NEWTON. FOR. 
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The  Histogram  Method  for  the  Solution  of  Welbixll  Parfl"»«»ter8 
Three-parameter  Welbull  distribution  is  expressed  as: 


P(0) 


1  -  exp  I  -c(a  -  Oy) 


} 


(c.l) 


The  function  F(0)  can  be  determined  in  several  ways.  In  the  mean 
raiik  method  n  tests  «ire  conducted  and  the  data  is  ordered  by  increasing 

fracture  stress.  This  yields  a  sequence  of  numbers  for  Oj  < . <  a^<  a 

The  estimation  of  cumulative  distribution  function,  F(a) ,  is  then  given 
by 

F(a.)  =  - -  ;  also  FCa.)  =  f^  (C.2) 

^  n  +  1 


A  more  sophisticated  statistical  treatment  [38]  known  as  the  Median 
Rank  method  uses  the  estimator 


i  -  .03 
n  +  .04 


(C.3) 


This  formula  is  an  approximation  to  Median  Rank  values  from  the 
imcomplete  Beta  function. 

Other  techniques  exist  vrtiich  also  seek  to  convert  a  set  of  data, 
i.e.,  a  number  of  specimens  with  a  number  of  fraction  values  to  a 
probability  of  fracture  table.  Basically  this  involves  the  integration 
of  a  histogram.  Figure  C.l, or  frequency  distribution  P(o),to  obtain  the 


Stress  Level 


Stress  Level 


Figure  C.l.  Probability  Distribution 


Thus  P(in  the  range  0^  to  ) 


0^  +  The  Number  of  failures  in  the  range  0^  to  o^+i 


And  once  P(a)  is  known. 


0 

F(a)  =  j  P(0)  do 

—  QC 


O 

j  p(a)  da 


(C.4) 


Numerous  techniques  exist  to  convert  this  histogram  information  to  a 
probability  of  failure  curve,  it  is  assumed  then  that  F(a)  will  be 
known  for  a  variety  of  data  points,  0  so  that  a  table  may  be  devised 
for  0  versus  F  .  Thus,  it  is  required  that  values  of  the  distribution 
constants  be  determined  for  the  F(o)  curve. 


A  typical  example  below,  T2d)le  C-1,  shows  how  to  obtain  a 


No.  of  Failures 


failure  probability  curve  from  the  discrete  element  (histogram) 
approach . 


Stress 

®i 

Rank 

i 

Probability  of  Failure  F(a) 

F(a^)  vs 

10.02 

1 

0 

0 

0  ^  a  <  10.02 

12.30 

2 

1/12 

0.083 

10.02  ^  a  <12.30 

14.76 

3 

2/12 

0.166 

12.30  <  14.76 

15.30 

4 

3/12 

0.250 

14.76  ^  a  <15.30 

16.24 

5 

4/12 

0.333 

15.30  ^  C  <  16.24 

16.94 

6 

5/12 

0.416 

16.24  ^  0  <  16.94 

17.20 

7 

6/12 

0.500 

16.94  <  17.20 

17.82 

8 

7/12 

0.583 

17.20  <  17.82 

18.40 

8/12 

0.666 

17.82  ^  a  <18.40 

19.60 

9/12 

0.750 

18.40  ^  a  c  10.60 

20.15 

10/12 

0.833 

19.60  :<a  <20.15 

22.40 

!■ 

11/12 

0.916 

20.15  iS'a  <22.40 

Table  C-1  .  Stress  I^vel  vs  Probability  .  "  Failure 


The  histogram  2md  the  cumulative  distribution  function  of  the  data 
presented  in  Table  C-1  is  displayed  as  follows : 


Figure  C2.  Cumulative  Frequency  Distribution  of  Table  C.l  data 


[•TiTi] 


Stress  Level,  a 


Figure  Q3.  Probability  of  Failure  for  the  data  of  Table  C-1. 

According  to  least-square  analysis  the  residue  equation  expressing 
the  sujn  of  the  square  of  differences  between  data  and  predictive 
equation  must  be  minimized.  The  residue  equation  with  weights  where 
the  cumulative  integrated  distribution  is  assumed, Eq. (C .4)  is  given 
by: 

22.>f0 

Res(c,m,Gy)  =  J  w(o)  jF(a)  -  F(C)  do  (C.5) 

10.02 


where  w(C)  =  weight  per  unit  stress. 


Expansion  of  Eq.  (3.5)  takes  the  form; 

/ 


Res(c,m,ay) 


12.jp 


J  w(a)  F(o)  do  + 


10.02 
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14  .  7E 

J  w(a){^-  F(a)}^  do  + 

12.30 

or  in  summation  form: 


n-1  __  2 

Res(c,m,a  )  =  Z  V{a^)  |f(Oj(^)  -  F(O^)  |  da  (C.7) 

i=l 


where  w(a)  is  assumed  cons tar t  in  each  range  expressed  in  Table  C-1. 

The  solution  of  Eq.  (C .7)  determines  the  values  of  three  Weibull 
parameters  c,  m,  and  However  this  method  is  veiry  cumbersome  and 

sometimes  the  predictive  equation  for  F(a)  is  not  of  closed  form  as 
is  the  case  of  4-point  bending  hollow  tube. 
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APPENDIX  D 

Statistical  Properties  of  Welbull  Distribution 


The  three-parameter  Velbull  distribution  In  Its  cumulative  form 
Is  given  as 


(D.l) 


For  the  case  =  0,  the  Eq.  (4.1)  takes  the  form  of  two-parameter 
Welbull  distribution  given  as: 


F(CJ) 


(D.2) 


In  this  case  the  fracture  can  take  place  at  any  positive  value  of 


the  fracture  stress. 

Since  the  three-parameter  distribution  can  always  be  converted 
to  the  two-parameter  distribution  by  a  single  linear  transformation* 
the  two-parameter  Weibull  distribution  is  used  to  illustrate  the 
properties  of  the  Weibull  distribution. 

The  probability  density  function  (p.d.f.)  is  obtained  by  differ¬ 
entiating  the  Eq.  CD.2)  . 


P(a)  = 


(D.3) 


taking  scale  par£uneter,  0=1,  the  Eq.  (4.3)  can  be  written  as: 

o 


m 

P(a)  =  jn(0)  exp(-(a)) 


(D.4) 


A  graphical  plot  of  the  Eq.  (D.4)  for  different  values  of 
shape  parameter,  m,  to  show  the  various  forms  of  the  probability 


density  function  has  been  shown  in  Figure  D-1. 

The  scale  parameter  Is  used  to  locate  the  Welbull  distribution 
along  the  O  eixis.  Assuming  a  >■  Into  the  cumulative  distribution 
function,  Eq.  (D.2),  It  can  be  easily  shown  that  for  any  Welbull 
distribution  the  probability  of  failure  prior  to  Cf^  Is  equal  to  63.2% 
cuid  Independence  of  the  shape  parameter,  m. 

F(a«a^)  =  1  -  exp 

m 

=  1  -  exp(-l) 


=  0.632 

From  the  above  fact  It  Is  quite  clear  that  the  scale  parameter, 
0^  will  always  divide  the  area  under  the  p.d.f.  Into  63.2%  for  all 
values  of  m. 


Figure  0.1.  The  Welbull  Distribution  for  Different  Values  of  m 


A  Goodness-of-Flt  Test  for  the  Welbull  Distribution 

When  dealing  with  Weibull  subjected  experinental  data  one  often 
faces  a  question  as  to  whether  the  data  fits  the  two-parameter  or 
three-parameter  Weibull  distribution.  The  sinplest  technique  one  can 
think  of,  is  to  transfer  the  tvns-parameter  Weibull  distribution  into 
a  linear  relationship  by  taking  the  log  twice  and  plotting  the  data 
on  a  graph  paper.  If  the  resulting  shape  of  Weibull  plot  is  a  straight 
line  then  one  can  eeisily  decide  that  the  data  follows  a  two-parameter 
Weibull  distribution.  Failure  to  get  a  straight  line  plot,  and 
especially  when  the  plot  curves  downwfurds  at  the  lower  end,  gives 
the  indication  that  the  experimental  data  follow  a  three-parameter 
Weibull  distribution.  In  this  case  the  value  of  the  threshold  stress 
Oy  must  be  calculated  to  get  a  straight  line  plot. 

Mann  et.al.  [39^  developed  a  goodness-of-fit  test  especially  for 
the  Weibull  distribution  which  is  expected  to  be  more  powerful  than 
any  of  the  general  goodness-of-fit  tests. 

The  null  hypothesis  in  this  case  is  that  the  experimental  data 
is  two-p^u:ameter  Weibull  distribution.  If  the  null  hypothesis  is 
rejected,  then  other  distributions  including  the  three-parameter 
Weibull  distribution,  should  be  considered.  The  mathematics  of  the 
goodness-of-fit  test  is  quite  single  and  is  given  as  follows: 

Let  0^,0 . ,C^  represents  the  first  r  ordered  failure 

stresses  resulting  from  placing  n  specimen  on  test  and  truncating  the 
test  at  the  time  of  r^  failure  (r^n) . 


,r.,  the  test 


De£inin9  EK^  as  *  InCo^)  for  1*1*2/ . 
statistics  is  given  as; 


r-1 

£ 

i=(r/2)+l 


Mi 


r-1 

Z 

iwl 


(EKj_^.i  -  EK^) 


Mi 


(D.5) 


where  (r/2)  denotes  the  greatest  integer  r/2j  for  example,  if  r  *  9 
then  {r/2)  ®  4.  The  value  for  Mj|^*s  are  found  in  l38),  along  with  the 

critical  value  for  S . 


UNI AX. FOR 


C  THIS  PROGRAM  IS  FOR  THE  UNIAXIAL  CONSTANT  STRESS  FIELDS 
C  PROGRAM  WEIBUL  THREE  PARA.METERS 
C*  THE  PARAMETERS  ARE  C,  M,  AND  SIGU 
C  THIS  PROGRAM  COMPUTES  ALL  PARAMETERS  BY 
C  LINEAR  REGRESSION  METHOD.  FOR  NOTATION  SIGU-SIGO 
C  THIS  PROGRAM  HAS  BEEN  EXECUTED  ON  DEC-1060  SYSTEM 

DIMENSION  X(20),F(2C),W(20),EK(20),P(20),WN(20),WM(20) 
WRITE{5,10) 

10  FORMATC IX, "ENTER  DATA  FILE  NO.*) 

READ(5,*)  ND 
READ(ND,*)  N 

C  ND=NUMBER  OF  DATA  FILE 
C  N=NUMBER  OF  WEIBUL  POINT  PAIRS 
C  W-ARBITRARY  WEIGHTS 

X-FRACTURE  STRESS  VALUES 

READ(ND,*)  ((X(I),F(I),W(1)),1=1,N) 

DO  1  1=1, N 

1  WR1TE(5,20)  (1,X(I),F(I).W(I)) 

20  F0RMAT(1X,*X(I),F(I).W(I)  FORI=* , 12, 3(1X,E12.6) > 

DO  2  1=1, N 

EK(I)=AL0G(AL0G(1./(1.-F(1)))) 

2  CONTINUE 
NN=5 
JK=0 
1NT=10 
RESO-l.E+30 
BEG=-.99*X(1) 

EMD=.99*X(1) 

98  H=(EMD-BEG)/FLOAT(INT) 

DO  95  1=0, INT 
SIK=BEG+FLOAT(I)*H 
DO  96  11=1, N 

96  P(ll)=ALOG(X(II)-SIK) 

CALL  LIN(W,WN,P,EK,N,S1K,A,B,RES) 

IF(RES.GT.RESO)  GO  TO  95 

RESO=RES 

SIGU=S1K 

AG=A 

BG=B 

J»I 

95  CONTINUE 

WRITE(5,80)  AG,BG,RESO,S1GU,SIGO 
C  REDIFINE  EMD  AND  BEG  AND  CONTINUE 
IF  (JK.GT.NN)  GO  TO  97 
JK-JK+1 

IF  (J.EQ.O)  J=1 

IF  (J.EQ.INT)  J-INT-1 

EMD=BEG+FLOAT( J+1 )*H 

BEG=EMD-2.*H 

RESO=RES 

SI GO- SIGU 

GO  TO  98 


80 


yu 


FORMAT(2X,*AG»*,E12.6,7X,*BG«*,El2.6/,2X,‘’Ri:SO=>",El2.6, 

1  5X.’SIG'J“*,E12.6,5X,*SIG0»‘‘,E12.6) 

C  EC.ESIGO,  AND  EM  ARE  THREE  WEIBULL  PARAMETERS 
97  EC»EXP(AG) 

ESIGO-SIGU 

EM-BG 

VOL»l 

C*  VOL  SHOULD  BE  CHANGED  AS  PER  PROBLEM 
VOLUN-1 

SIGO=>(1./(VOLUN*EC))**(1/EM) 

WRITE(5,100)EC,EM,SIGU,S1G0 
100  FORMAT (/,20X.‘’WEIBUL  PARAMETERS  ARE“,//,2X, 

1  •EC=*'’,E12.6.5X,‘'EM=‘',F8.2,5X,'’S1GU='’,E12.6,5X,*SIG0='’,E12.5,/) 

C  CALCULATED  PROBABILITIES  ARE  LISTED  AS  BELOW 
DO  130  11-1, N 

130  P(II)“AL0G(X<I1)-SIGU) 

CALL  LIN(W.WM.P,EK,N,SIGU,A,B,RES) 

DO  108  1=1, N 

PCAL=1.-EX?(-(VOL/VOLUN)*((X(I)-ESIGO)/SIGO)**EM) 

108  WR1TE(5,110)  PCAL.F(I),W(I) 

no  F0RMAT(1X,‘‘PCAL.F(I),W(I)='*,E12.6,2(1X,E12.8)) 

STOP 

END 

C  SUBROUTINE  LIN  CALCULATES  VALUES  OF  CONSTANTS  A  AND  B 
SUBROUTINE  LIN( W, WN, ?, EX, N, SIGU, A, B, RES) 

DIMENSION  W(20),P(20),EK(20),WN(20) 

C«0. 

D»0. 

E=0. 

G-0, 

H=0. 

DO  4  I"1,N 
C=C+(W(I)) 

D'=D+(W(I)*P(1)) 

E-E+(W(I)*?(I)*P(I)) 

G'=G+(W(I)*EK(I)) 

H»H-t-(W(I)*EK(l)*P(l)) 

4  CONTINUE 
DEN=E*C-D*D 
A»((E*G-D*H)/(DEN)) 

B»((C*H-D*G)/(DEN)) 

C  THE  CONSTANTS  A  AND  B  ARE  KNOWN 
C  NOW  FIND  RESIDUE 
SUMM=0. 

DO  5  1-1, N 

WN(I)"(A+B*?(I))-EK{I) 

5  SUMM“SUMM-*-W(l)*WN(l)*WN(l) 

RES--*SUMM 

RETURN 

END 


EX  UNIAX.FOR 

FORTRAN;  UNIAX  Ql 

MAIN . 

LIN 

LINK:  Loading 

[LNXXCT  UNIAX  Execution] 

ENTER  DATA  FILE  NO. 

X(I),F(I).W(I)  FORI=  1  .755000E+01  .625G00E-01  .lOOOOCE+Gl 

X(I).F(I),W(I)  FORI-  2  .855000E+01  .125000E+00  .lOOOOOE+01 

X(I),F(I),W(I)  FORI-  3  .109100E+02  .137500E+00  .lOOOOOE+01 

X(I).F(I).W(I)  FORI-  U  .11600CE+02  .250000E+00  ,  .lOOOOOE+01 

X(I).F(I).V.'(I)  FORI-  5  .UAiOCE-^O’  .31250CE+00  .lOOOOOE+01 

X(I),F(I),W(I)  FORI-  6  .135000E+02  .375000E+00  .lOOOCOE+01 

X(I),F(I),W(I)  FORI-  7  .136500E+02  .A3750CE+00  .lOOOOOE+01 

X(I),F(l),wa)  FORI-  8  .140000E+02  .5CC000E+00  .lOOOOOE+01 

X(I),F(I).W(I)  FORI-  9  .1A7500E+02  .562500E+00  .lOOOOOE+02 

X(I),F(1).U(1)  FORI-10  .148300E+G2  .625C00E+00  .lOOOOOE+01 

X(I).F(I).W(1)  FORI-11  .159600E+02  .6S7500E+00  .lOGOOOE+01 

X(1),F(I).W(I)  FORI-12  .16910CE+02  .75000CE+00  .lOGGOOE+01 

X(I),F(1),W(I)  FORI-13  .16600GE+02  .8125CCE+00  .lOOOOOE+01 

X(I),F(I).W(I)  FORI-14  .2QOOO0E+O2  .875000E*00  .lOOOOOE+01 

X(I).F(I),W(I)  FORI-15  .215000E+G2  .932500E+00  . lOOOOOE+01 

AG--.103722E+02  BG-  .377017E+01 

RESO-  .682176E+00  SIGU-  .OCOOOOE+00  SIGO-  .OOOOOOE+00 

AG— .103722E+02  BG-  .377C17£+0i 

RESO-  .632176E+00  SIGU-  .^A7035E-07  SIGO-  .OOOOOOE+00 

AG--.103722E+02  BG-  .377017E+01 

RESO-  ..‘I82176E+00  SIGU-  .558794E-07  SIGO-  .447G35E-07 

AG— .103i60E+02  BG-  .376275E+01 

RESO-  .88217i*E+00  SIGU-  .239185E-01  SIGO-  .538794E-07 

AG--.103i60E+02  BG-  .376275E+01 

RESO-  .682174E+00  SIGU-  .239135E-01  SIGO-  .239185E-01 

AG--.103i^GE+02  BG-  .376217E+01 

RESO-  .682173E+00  SIGU-  .258319E-01  SIGO-  .239185E-01 

AG--.103436E+02  BC--  .376207E+01 

RESO-  .S82173E+00  SIGU-  .261189E-01  SIGO-  .258319E-01 

WEIBUL  PARAMETERS  ARE 

EC-  .321973E-04  EM-  3.76  SIGU-  .261189E-01  SIGO-  0.15634E+02 

PCAL,F(I),W(1)-  .618403E-01  .625000E-01  .lOOOOCE+01 

PCAL.F(I),W(I)-  .97G441E-01  .125000E+00  .lOOOOOE+01 

PCAL,F(I),W(I)-  .225875E+0C  .137500E+00  .lOOOOOE+01 

PCAL,F(I),W(I)=  .275760E+00  .250000E+00  .lOOOOOE+01 

PCAL,F(I),W(1)-  .342910E+00  .312500E+00  .lOOOOOE+01 

PCAL,F(I),W(I)-  .435352E+00  .375000E-00  .lOOOOOE+01 

PCAL,F(1),W(I)-  .448912E+00  .437500E+00  .lOOOOOE+01 

PCAL,F(I),W(I)=  .480830E+00  .500000E+00  .lOOOOOE+Ol 

PCAL,F(I),W(I)-  .549775E+00  .562500E+00  .lOOOOOE+02 

PCAL.Fd) ,W(I)=  .557114E+00  .625000E+00  .lOOOOOE+01 

PCAL.Fd)  ,W(I)=  .658394E+00  .68750GE+00  .lOOOOOE+Ol 

PCAL,Fd),Wd)=»  .696612E+00  .75C0O0E+O0  .lOOOOOE+01 

PCAL.Fd)  ,Wd)-  .712237E+00  .812500E+00  .lOOOOOE+01 

PCAL.Fd), Wd)=  .919005E+00  .875000E+00  .lOGOOOE+01 

PCAL.Fd)  ,Wd)-  .963132E+00  .932500E+00  .lGOOOOE+01 

STOP 

END  OF  EXECUTION 

CPU  TIME:  1.69  ELAPSED  TIME:  2:13.53 

rvTT 
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BENDIN.FOR 

C  THIS  PROLRAM  IS  FOR  THE  4-POlN’T  LOADING  PURE  BENDING  STRESS 
C  PROGRAM  WSIBUL  THREE  PARAMETERS 
C  THE  PARAMETERS  ARE  SIGMAO,  M  .AND  C 
C  THIS  PROGR.VM  COMPUTES  ALL  PARAMETERS  BV 
C  LINEAR  REGRESSION  METH'^D.  FOR  NOTATION  SIGMA0»SIG0 
C  THIS  PROGRAM  HAS  BEEN  EXECUTED  ON  DEC-1060  SYSTEM 

DIMENSION  X(20).F(20).W(20),EK{20),P(20),WN(20),WM(20) 
wr:te(5.io) 

10  FORMAT( IX, "ENTER  DATA  FILE  NO.*) 

READ(5,*)  ND 
READ(ND,*)  N 

C  ND=NUMBER  OF  DATA  FILE 
C  N-NUMBER  OF  WEIBUL  POINT  PAIRS 
C  W=ARBITRARY  WEIGHTS 
C  X=FRACTURE  STRESS  VALUES 

READ(ND,*)  ((X(I),F(I),W(I)),I=1,N) 

DO  1  1=1. N 

1  WRITE(5.20)  (1.X(I),F(I),W(I)) 

20  F0RMAT(1X,°X(1),F(I),W(I)  F'^R1=" ,12,3(1X,E12.6)) 

DO  2  1=1. N 

EK(I)=AL0G(AL0G(1./(1.-F(1))))+AL0G(X(I)) 

2  CONTINUE 

C  V  IS  THE  VOLLE-IE  OF  TEST  SPECIMEN 

V=1.0  ! SHOULD  BE  CHANGED  AS  PER  VOLUME 

NN=5 

JK=0 

INT=iO 

RESO=1.E+30 

BEG=-.99*X(1) 

EMD=.99*X(1) 

98  H=(EMD-BEG)/FL0AT(1NT) 

DO  95  1=0. INT 
SIK=BEG+FLOAT(I)*H 
DO  96  11=1, N 

96  P(1I)=AL0G(X(1I)-SIK) 

CALL  LIN(W,WN,P,EK,N,SIX,A,B,RES) 

1>(RES.GT.RES0)  GO  TO  95 

RESO=RSS 

SIG=SIX 

AG=A 

BG=D 

J=I 

95  CONTINUE 

WRITE(5,80)  AG.BG.RESO.SIG.SIGO 
C  REDIFINE  EMD  AND  BEG  AND  CONTINUE 
IF  (JK.GT.NN)  GO  TO  97 
JK=JK-t-l 

IF  (J.EQ.O)  J=1 

IF  (J.EQ.INT)  J=INT-1 

EMn"BEG+FLOAT(J+l)*H 

BEG=EMD-2.*H 

RESO=RES 
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SIGO-SIG 
GO  TO  98 

80  FORMAT(2X,‘‘AG=“,El2.<i,7X,*BG»"  ,E12.6/,2X,*RESO»“  ,E12.6, 

1  5X,*SIG»‘*,E12.6,5X,*SIG0«',E12.6) 

C  EC.ESIGO,  AND  EM  ARE  THREE  WEIBULL  PARAMETERS 
97  SIGG— ( (AG-AL0G(V/2)+AL0G(BG) )/( BG-1 . ) ) 

SIGO=>EXP(SIGG) 

ESIG0=S1G 
EM^BG-l . 

WRITE(5.100)SIG0,EM,S1G 

100  FORMAT(/,2X/ WEIBULL  PARAMETERS  ARE :  ‘  ,2X,’S1GMA0«‘’ ,E12 .6 , 

1  5X.‘’iii=-‘’,El2.5,5X,‘’SIGU=‘’,E12.6/) 

C  CALCULATED  PROBABILITIES  ARE  LISTED  AS  BELOW 
DO  130  1I«1,N 

130  P(1I)=AL0G(X(1I)-S1G) 

CALL  LIN(W,WM,P,EK,N.SIG,A,B,RES) 

DO  108  1=1. N 

PCAL=1 .0-EXP(-(V/(2*(EM+l .)))*( 1 .-(ESIGO/X(l) ) ) 

1  *((X(I)-ES1GO)/S1GO)**EM) 

108  WRITE(5,110)  PCAL.F(1),W(I),WM(I) 

110  F0RMAT(1X.“PCAL,F(1),W(I),WM(1)=‘’,E12.6,3(1X,E12.6)) 

STOP 

END 

C  SUBROUTINE  LIN  CALCULATES  VALUES  OF  CONSTANTS  A  AND  B  WHICH 
C  REPRESENTS  A  STRAIGHT  LINE  V=A+BX 

SUBROUTINE  L1N(W,WN,P,EK,N,SIG,A,B.RER) 

DIMENSION  W(20) ,?(20) ,EK{20) ,WN(20) 

C=0. 

D=0. 

E=0. 

G=0. 

H=0. 

DO  4  1=1, N 
C=C+(W(1)) 

D=D+(W(I)*P(I)) 

E=E+(W(I)*?(I)*P(I)) 

G»G+(W(I)*EK(I)) 

H=H+(W(I)*EK(I)*P(I)) 
k  CONTINUE 

DEN=E*C-D*D 

A=((E*G-D*H)/(DEN)) 

B=((C*H-D*G)/(DEN)) 

C  THE  CONSTANTS  A  AND  ^ARE  KNOWN 
C  NOW  FIND  RESIDUE 
SUM!-1=0 . 

DO  5  I-1,N 

WN(l)»(A+B*P(i: )  )-EK(  I) 

5  SUMM=SUMM+W(I)*WN(1)*WN(I) 

R£S=SUMM  ^ 

RETURN 

END 
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EX  BENDIN.FOR 
FORTRAN;  BENDIN 
MAIN. 

LIN 

LINK;  Loading 
[LNKXCT  BENDIN  ExacucionJ 

ENTER  DATA  FILE  NO. 

39 

X(I),F(I),W(I)  F0RI=  1  .lOGGGOE+02  .3005SOE+00  .lCOOOOE+01 

X(I),F(1),W(I)  FORI=  2  .120G0GE+02  .469116E+00  .lOOOOOE+01 

X(I).F(1),W(I)  FORI=  3  .15GGGGE+02  .6990CCE+00  .lOOCOOE+01 

X(I),F(I),W(I)  FORI=  4  .170CCGE+02  .814079E+00  .lOOOOOE+01 

X(I),F(I),W(I)  FORI=  5  .20000GE+02  .92272GE+00  .lOOGOOE+01 


AG=-.555985E+01 

BC-= 

.328S26E+01 

RESO=  .54Q424E-03 

SIG-- 

.1980GCE+01 

SIGO=  .OOOCOCE+00 

AG=-.440115E+01 

BG= 

.295214E+01 

RESO=  .182929E-04 

SIG= 

.316800E+G1 

SIGO=  .198G00E+01 

AG=-.455A4CE+01 

BG= 

.299723E+01 

RESO«  .614190E-07 

SIG= 

.300960E+G1 

SIGO=  .316800E+01 

AG=-. 4569745+01 

BG= 

.3C0175E*01 

RESO=  .321390E-G7 

SIG= 

.299376E+G1 

SIGO=»  .300960E+01 

AG=-.A5636GE+G1 

BG= 

.2999955^01 

RSSO=  .60^9575-03 

SIG= 

.3G0G10E+01 

SIGO=  .299376E+01 

AG^-.456360E+01 

BG= 

.299995E*01 

RESO’’  .604957E-08 

SIG= 

.30GG10E-^01 

SIG'^=  .300010E-*-01 

AG=-.45634S5+01 

BG*- 

.2999915+01 

RESO*  .602G20E-08 

SIG= 

.3G0G22E+G1 

SIGO=  .300010E+01 

WEI3ULL  PARAMETERS 
SIGMA0=*  .399656E+01 

ARE: 

tn= 

0.199995+01 

SIGU=  .300022E+01 

PCAL,F(I),W(I),WM(I)=  .300532E-00  .30G580E--00  .lOQOOOE+01  .891089E-05 

?CAL,F(I),W(I).WM(I)='  .469106E+a0  .469116E+G0  .ICGOOOE+Ol  -.289977E-C4 

PCAL,F(I),V(I),WM(I)=  .699021E+00  .G99000E+G0  .190C0CE+01  .574S87E-C4 

PCAL,F(I),'N(I),V--M(I)=  .SIAGGGE-^OO  .814079E+00  .lOOOOOE-^01  -.418636E-04 

PCAL,F(I),W(I),v;m(I)«’  .922722E^00  .922720E+00  .lOGOOOE+01  .768900E-05 

STOP 

END  OF  EXECUTION 

CPU  TIME;  0,94  ELAPSED  TIME:  1;5,62 
EXIT 


HOLLOW. FOR 


VROCR^X  BEI'.niK 

0*  THIS  RROCRAM  IS  FOR  A-FOliST  BEKniKC  HOLLOW  CERAMIC  ROI> 
REAL  L/lCAC) 

DI MEilS lOi:  ER F C  AC  )  , TFF (  AO )  >’( 40 )  ,  AA (  5  )  ,  y  ( 4C  ) 

COMMOW  / y AR / 1 , L , C , M , El , S ICO , EM , RO , RI 
WRITE (5, AC) 

AO  FORMATC^X/EIITER  RATA  FILE  HOOMBER  AIID  ERROR  BOUJiD") 

R,EAn(5,*)0n,KRR 
C*  K*=i:o.  OF  DATA  FOIilTS 
C*  y*LOAn  VALOES 

REAP ( on, ■^)  L,RO,RT 

READ(;;n,*)K 

KS>=] 

REAnCOD,*)  y(.i),w(i) 
no  12. 

RKAOCif),-^)  ?(1)>’(I) 

IF  (y(I).EQ.y(I-!))  CO  TO  22 
KS"=KS>1 

22  COOTIOUE 

X*''L0AT(KS<-1  ) 

C*  i.'CW  nC'  PACK  METHOD 
J*1 

r’yF(l)«l  ./X 

nu  t  i‘'2,K 

r;y!'(  1  )*=i:y‘-'(  i-i ) 

IF  (y(I).r;Q.'^(I-l ))  Gu  T<.  1 
J*  T  +  1 

r  RFC  I)-=F1.0ATf  J)/X 

5  co;:ri;:;jE 


;:(;w 
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OP  MAT  (1 X ,  “  J<=  ’  ,  1 2 ,  ’  E  PF  f  J ) «  ’  ,  7 . 5  ) 

W?.ITrri,Al  ) 

Fop matc  i  x  ,  ’  eotf ?  i. ,  p.g ,  ri  ,  err  ’ ) 

CALCULATE  MoMi  OT  FRi.-M  I*,‘A0  VALOES 
no  2  I*=1,K 
M(  1)>=V(1)-v'j.h/2. 

.  '•J./A.)+(Ro-'fVi-px*-*A) 

WP  !  T '  ( =5 ,  Q  )  1. ,  P.O ,  P I ,  ERR ,  r  I 

Fu^MAI(//^X,’L>=^1X,F7.?,?X,^?l>^1X,FS.3,'iX,»R>’,lX 

r  ^  ,  PX  ,  ■’  r  PP-'’ ,  1 X  ,  '’7 .  -7  ,  PX ,  "  El-'’ ,  1 X  ,  F7  .  A  ,  /  ) 

’.-ODr  ;:0.  Cr'TTlirP  1  i,?  2  OP  '' ) 

^CE  Lr.;:CT!l  OF  S'^-ClMr;; 

.;e-i  (c  va'’!;-')  ;;;;«=2  fEM  varies)  '.lo-p  Csico  vap 

WRITE  ('),•*  2  ) 

!'0P'-  '.Tn:<,'’i:;y;:T  c,EM,sic:',.;;:(conr )  ,ft;;al,;:os'’ ) 

AAf  !  )>=0 


AA(  2  )*=E  '' 

■.A(  - 

)  Ci,  lo 
.)G'  2v  -o 

atc.r  n 

Tt  ;;  '  j  : 
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DO  17  1*1, K 

vap,-=aa(n;:)+float{ij-i  )*i<:i;c 

N*4 

IF  (t;K.r:Q.l)  C*VAR 
IF  (;m.KQ.2)  EM*VAR 
IF  (i:n.EQ.3)  SICU«VA? 

CALL  ACAR(f:PR,iJ,BK,XX) 

TPF(I)*XX 
1 7  CONTIiJIJE 

C*  :;OW  COMW'JTE  residues 
RES«0. 

DO  5  1*1 ,K 
WU*(EFF(1 )-TVF(l)) 

5  RES*RrS+W(I)*WU*W;: 

1 F  ( Nil .  E Q .  1  )  WRITE  (  5 , 1 5  )  VAR ,  M ,  S I CU ,  P  r  S 
IFfr.U .  EQ.  2)WRITE(5 , 16)C,  VAR ,  SICU ,  RrAS 
IF  ( i:X .  EQ .  3  )  WRITE  (  5 , 1 8  )  C ,  EM ,  VAR ,  RE  S 
1 5  FORMAT  ( 1 X ,  “  V AP. :  C-^%  E 1 4 . 6 , 2X , EM*  “  ,  1 X ,  F  6 . 3 , 2X ,  “  S 1CU-:  “  , 

1  lX,ElC.3,2X,"RES*'’,lX,ElA.b) 

It)  F0RMAT(1X,’C*'’,1X,E13.7,1X,'’VAR:EM*‘*  ,1X,FS.3,1X/S1G:;*’, 

1  lX,KK’.'i,lX,''?.ES*'*,lX,E].*.6) 

1 3  FORMATC  IX , '=C*“ ,  1 X ,  El  2 . 7 ,  IX ,  ’ EM*'’ ,  IX , F6 . 3 ,  IX ,  “ VAR:  SIC:J*“ ,  1 X, 

1  E14.fc.,lX,’?ES**,lX,rl4.6) 

IF(  IIOS  .UK . C)  CO  T(j  )' 

no  81  1*1 ,K 

WRITE(5 , 7<^ )  ( 1  ,EPF( I  )  ,TyF(  1 )  ,W( T  ) ) 

70  FOPMATClX,’ryFCr),TPF([;,W(’T)  FOR  I*’ ,  12, 3(2X,  E 1 3 . 3)) 

81  uOATIixiE 

IL  CCUTIUUE 

CO  TO  1 1 
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SUBF.O'o'TIIlH  ACARCiiRP  ,K,1U;,XXX) 

RKAL  L,M(u(;) 

COMMON  /PAP/  l,L,C,M,F.T,SIC!j,r:M,PO,RI 
C*  SIMPSON  RULK  PROGRAM  FOR  ClRCUl.AR  BKNniKC  GASF 
C*  A  ANB  B  ARE  LOVER  AND  UPPER  LIMIT  OF  INTECPATTOf. 

C*  FOR  FIRST  PART  OF  INTECPATION 

C*  D  AND  A  AP,E  LOWER  AND  UPPER  LIMIT  OF  INTTECRATION  FOR 
C*  SECOND  FART  OF  INTEGRATION 
C*  N  IS  THE  NUMBER  OF  SUBDIVISION 
C*  DX  IS  THE  INTERVAL 

C*  C, EM, AND  SIGN  ARE  MATERIAL  PARAMETERS 
FN-=FLOAT(N) 

BNT-=0 . 

C+  EI^MO'LE]:!  OF  INERTIA 
C*  M*BEr:DINC  MOMENT 
C*  D  IS  y-THRESHoLD 

C*  TIIPEE  CASES  HAVE  BEEN  CONSIDERED  FOR  Y-THP.ESHOLD 
A*RI 


B«=Ro 

D*SICU*EI/M(I) 

C*  CASE-1.  Y-THRESHOLD  CAN  OCCUR  (MJTSTDE  PO.CIF  OCCURS, THE!; 
C*  CASE-2.  Y-THRE'^HOLD  CAN  OCCUR  BETWEEN  O.AND  PI 
C*  CASE-3.  Y-THRKSHOLD  CAN  OCCUR  BETWEEN  RI  AND  RO 
] :  0  I F  ( D .  C  E .  U  .  A  N  D .  D .  LT .  A  )  C  A  LL  S I M  V  ( A ,  B ,  D ,  F  N ,  B  N  ) 

I F  ( D .  C  T .  A .  AN  D .  D .  LT .  B  )  C  A  LL  S 1  MM  ( A ,  B ,  D ,  K  N ,  B  N ) 
IF(n.CT.B)  CO  TO  3 

C*  XXX^nirRtjTlVAl  PROBABILITY  oF  FRACTURE 
XXX*=l.-rXP(-Bi;) 

iFf ab'=(b’;-b::t).le.err)  return 
fN'=f::*2 

oO  To’  1  W'' 

3  BN“=;  . 

XXX«='..  . 

Rr.Ti'R.T 
END 


l,l.,C,M,rI,SICU,EM,Ru,RI 


SUBROUTINE  SlMy(A,B,D,F:;,BN) 
REAL 

REAL  L,M(4;.) 

COMMO"  /PA?/ 

Il^IFIXC^N) 
nX’*:(B-A)/FN. 
dxu^Ca-o/f;; 

T1'.X1*’’<’.-‘NLX1 

r::x2’^:'  .*v-x-: 

rNll-GION;’’ !  fAl^CRAN'UH) 
r:;:l2’'CRANFUfN)-’-CRAN’-2(A) 
lALlS:'  r.NN  ’  NN  E::3 


;)i  1-1 

Yl'-..*-:  Xl-’-U  ATM) 
Y.:-r+NX2 "'U’  './(j) 


,  BN-=0 ) 


i 


5 


'JOIITIKUK 

Y]*A 

Y2*D 

DO  4  J*J  ,iJM 
Yl-YUTDXl 
Y2-Y2+TDX2 
i!:i;31*KN3]+GRAi;FJ(Y]) 

F:::32*r::;32+CRAKF2(Y2) 

Bt5«(nXl*(E:;ll+4.*£N21+2.*F.!333  )/3  . )  +  (DX2*(ED]  2+4 .*£1322 
+2.*r::)32)/3.) 

RKTURN 

EiiD 

I'UliCTIOi;  CRAKFUY) 

RrAl.  L,M(4C) 

COMMO!.'  /VAR/  l,l.,C,M,ET,SICIJ,r:M,PO,RT 
SIC*=M(I)*Y/£I 

I»'((f;iC-FlC3J).l.T.0)RTG«=SICT] 

rORMAT(2X/M,Y,SIGIJ,SIC,fVHn!;  STG-SIGIJ.LT.CO’''’ ,  ^X,4r  11.4) 
CRAi3F]‘'(2.*L*t;)*(((M(I)*Y/fc:i)-SIGi;)**£M)* 

( (RO+RO-Y+y)**.^) 

RKTUPi; 


EiiD 

FUGOTio;:  cra;:f2(y) 

REAL  l.;'‘U4A) 

COMMO;:  /RAR/  T,L/:,M,£I,SICU,r.M,RO,P.l 
SIC»=Mfl)*Y/r.T 

IF  ( (  S.IC-'^  ICIJ)  ,I.T .  I  )«1G«='^1C;J 

GTI  Fo?MAT(2X,''M,Y,SIG:',SIG,(VHi':'  SlORICr.lT.I'  i:3  C’^2)«'’,4r:] 

GPA:;I'2*=(2*L-^C)-*^(  f  (MCD+Y/KD-Sin.O-^+rlM)* 

3  ( ( ( ?0*RU-Y*Y.V** .  ^  )-f  ( RI*RI-Y*Y)** .  S )  ) 

PrTIjR;: 

? : JB  ?  o;;t  1  ■ :  r  i  mm  (  a  ,  b  ,  d  ,  i : ,  b  ;: ) 

REAL 

r-xi-=( 

Tr,xi>';’.+i)xi 

rxii  ■=c?.A:;r.i  (rO+GRAFFi  (  b) 


K;;2!  . 

.'■■iii'i  'I'. 

:x'  A 

Yl'=D*nX!+'^I,!'AT(J) 
f.;:2i*-'';;2i+GRA;:ri(Yi ) 


Yi-r 

Taj  4 

Yi'=Y’.  +t;-.xi 

r;:  !  (yi  ) 

\  \  *2 )/3.) 


M 
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NEWTON. FOR 


PROGRAM  NEWTON 
C*  NEWT0N-RAPH30N  METHOD 
C*  PROGRAM  WEIB’JL  THREE  PARAMETERS 
C*  THE  PARAMETERS  ARE  C,  M,  AND  SIGU 
C*  THIS  PROGRAM  COMPUTES  ALL  PARAMETERS  BY 
C*  NEWTON-RAPHSON  METHOD 

DIMENSION  X(20).F(20),W(20),A(20,20),Z(20).B(20) 
DIMENSION  U(20),U1J(20,3),U1(20),U3(20),UIKJ(20,3,3) 
WRITE(5,10) 

10  FORMAK IX, ‘ENTER  DATA  FILE  NO.*) 

READ(5,*)ND 

READ(ND,*)N 

C*  ND*=NO.  OF  DATA  FILE 

C*  N=  NO.  OF  WEIBUL  POINT  PAIRS 

C*  W=  ARBITRARY  WEIGHTS  ASSIGNED  TO  DATA  POINTS 

C*  X=  FRACTURE  STRESS  VALUES 

READ(ND,*)((X(I),F(I>,W(I)),I=1,N) 

DO  1  1=1, N 

1  WRITE(5,20)(I,X(I),F(I),W(I)) 

20  F0R^UT(1X,‘■X(1),F(1),W(1)  FORI=‘‘ ,I2,3(lX,E12.6)) 

C*  NOW  SET  THE  VARIOUS  SUMMATIONS  FOR  THE  LEAST  SQUARE 
C*  METHOD,  BUT  FIRST  MAKE  A  GUESS  FOR  C,  SIGU,  AND  M 
12  WRITE(5,24) 

24  F0R:UT( IX, ‘INPUT  THE  GUESSES  FOR  C,SI.GU,EM‘) 

READ(5,*)C,SIGU,ZM 

11  DO  3  1=1,  N 
T=(X(I)-SIGU) 

1F(T.LT.0.)  T=0. 

U(I)=G*T**EM 

S1=EXP(-U(I)) 

S2=S1*S1 

U1(I)=W(1)*((1.-F(1))*S1-S2) 

U3(I)=W(I)*(2.*S2-(1.-F(I))*S1) 

U1J(1,1)=U(I)/C 

UIJ{I,2)=-U(1)*EM/T 

UIJ(I,3)=U(I)*AL0G(T) 

UIKJ(I,1,1)=0. 

UIKJ(I,1,2)=U1J(I,2)/C 
UIXJ(I,2,1)=UIKJ(I,1,2) 
Ui:<J(I,2,2)=-(EM-l.)*UIJ(I,2)/T 
UIKJ(I,1,3)=UIJ(I,3)/C 
UIKJ(I,3,1)=UIKJ(I,1,3) 
UIKJ(l,2,3)=(l.+EM*ALOG(T))*UIJ(l,2>/EM 
UIKJ(I,3,2)=UIKJ(I,2,3) 
UIKJ(I,3.3)=U(I)*(AL0G(T))**2 
3  CONTINUE 

C*  NOW  FORM  THE  MATRIX  B(l),  AND  A(I,J) 

C*  ZERO  OUT  A  AND  B 
DO  30  1=1,3 
B(I)«0. 

DO  30  J=l,3 
A(I,J)=0. 
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1 

I  30  CONTIMUE 

DO  A  I»1,K 
DO  i  J=»1.3 

B(J)=»B(J)+U1(I)*U1J(I,J) 

DO  4  K»l,3 

A(J,K)“A(J,K)+U1(I)*UIKJ(I,K,J)+U3(1)*UIJ(I,J)*UIJ(I,K) 

4  CONTINUE 

CALL  SOLVE(A,Z.B.3,DET,5) 

WRITE(5,13)(Z(I).I=1,3) 

13  FOR.MAT(lX.“DC,DSIC-U,DM='‘’,E15.6) 

C=*C-Z(1 

S1GU=SICU-Z(2) 

EM=>EM-Z(3) 

WRITE(5,14)C.SIGU,EM 

14  FORMATCIX/NEW  WEIBULL  PARAMETERS  ARE:  C, S1GU,EM*=°  ,3E15 . 6) 

ITES1=1HY 

WRITE(5.15) 

15  FORMATdX/DO  YOU  WISH  TO  CONTINUE  WITH  PRESENT  CORRECTIONS®, 

1  /,® ANSWER  Y  OR  N“ ) 

READ(5,16)ITE 

16  for:-iat(ai) 

IF(ITE.EQ.ITESI)  GO  TO  11 
WRITE (5, 17) 

17  F0RMAT(1X,®D0  YOU  WANT  TO  STOP-ANSWER  Y  OR  N*) 

READ(5,16)ITE 

ifCite.eq.itesDgo  to  21 
GO  TO  12 

21  WR1TE(5,100)C,EM,SIGU 

100  FOR:iaT(/,2X,®C»'',E12.6,5X,“EM=®,E12.6,5X,®SIGU»*,E12.6/) 

STOP 

END 

SUBROUTINE  SOL VE ( A , X , B , N , DE T , N P ) 

DIMENSION  B(20),X(20),A(20,20),K(20),Y{20) 

DO  16  I^l.N 
16  K(I)-1 

N1=H-1 
DET=1. 

SIG=1 . 

DO  8  L“1,N1 

C  ****  SEARCH  FOR  LARGEST  ELEMENT 
C  NP=’UNIT  NUMBER  OF  PRINTER 

D»0. 

DO  1  L1=L,N 
DO  1  L2=L.N 

IF(ABS(A(L1,L2))-ABS(D))  1,1,1550 
1550  D»A(L1,L2) 

1D»=L1 

JD'’L2 

1  CONTINUE 
1F(D)2,99,2 

C  ****  INTERCHANGE  R^WS  AND  COLUMNS  TO  PUT  URGEST  ELEMENT  ON  DIAGONAL 

2  IF  (JD.EQ.L)  GO  TO  14 
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SIG—SIG 

1EMP»K(L) 

K(L)=*IC(JD) 

K(JD)=IEM? 

DO  U  I»1,N 
TEM?=’A(I,L) 

A(I,L)>»A(I,jn) 

A  A(I,JD)=TEMP 
lA  IF  (ID.EQ.L)  GO  TO  15 
SIG=-SIG 
DO  3  J=>L,N 
TEMP=>A(L,J) 

A(L,J)=A(ID,J) 

3  aUd,j)=temp 
TEMP=B(L) 

B(L)=B(ID) 

B(ID)=TEMP 
15  B(L)=B(L)/D 

C***EL1MINATE  ELEMENTS  IN  COLUMN  UNDER  LARGEST  ELEMENT 
L1=L+1 
DO  5  J=>L1,N 

5  A(L.J)=A(L.J)/D 
DET=DET*D 

DO  7  1=L1,N 
IF(A(I,L))  1515,7,1515 
1515  Di=Aa,L) 

DO  6  J=L1,N 

6  A(1,J)=A(I,J)-D1*A(L,J) 

B(I)=>B(I)-D1*B(L) 

7  CONTINUE 

8  CONTINUE 
1F(A(N,N))  9,99,9 

C***BACK  SUBSTITUTE  TO  SOLVE 

9  Y(N)=B(N)/A(N,N) 

DET=DET*A(N,N)*S1G 
DO  11  L»«1,N1 
LL«N-L+1 
D1=B(LL-1) 

DO  10  J=LL,N 

10  D1=D1-A(LL-1,J)*Y(J) 

11  Y(LL-1)=D1 
C****  RE-ORDER  ANSWER 

DO  12  I-1,N 
J^KCI) 

12  X(J)»Y(I) 

GO  TO  13 

99  WR1TE(NP,100) 

100  FORMAT(2X, 'MATRIX  IS  SINGUUR  NO  SOLUTI'^N  GIVEN') 
DET^O, 

13  RETURN 
END 
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DISTRIBUTION  LIST  FOR  H00019-T8-C-0520 


NO 


Commander 

Naval  Air  Systems  Command 
Washington,  DC  20361 
Attention:  AIR-OOD^  6 

AIR-320A  1 

AIR-5163D1+  k 

Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  2221? 

Attention:  Code  471 

Commander 

Naval  Surface  Weapons  Center 
White  Oak 

Silver  Spring,  MD  20910 
Attention:  Code  R31 

Director 

Naval  Research  Laboratory 
Washington,  DC  20375 

Attention:  Code  636O 

Commanding  Officer 

David  W.  Taylor  Naval  Ship  Research  &  Development  Center 
Annapolis,  MD  2l4l2 
Attention:  W.  Smith,  Code  2832 

Commander 

Air  Force  Wright  Aeronautical  Laboratories 
Wright-Patterson  Air  Force  Base 
Dayton,  OH  45433 

Attention;  Dr.  J.  Dill  POSL  1 

Dr.  H.  Graham  MLLM  1 

Brookhaven  National  Laboratory 

Upton,  NY  11973 

Attention;  Dr.  D.  Van  Rooyen 

Director 

Applied  Technology  Laboratory 

US  Army  Research  &  Technology  Laboratories 

Fort  Eustis,  VA  236o4 

Attention:  DAVOL-ATL-ATP  (Mr.  Pauze) 

US  Army  Research  Office 
Box  CM,  Duke  Station 
Durham,  NC  27706 

Attention:  CRCARD 


OF  COPIES 
11 


1 

1 

1 

1 

2 


1 

1 

1 
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Army  Materials  and  Mechanics  Research  Center 
Wai^ertown,  MA  02172 
Attention:  Dr.  R.  N.  Katz 

NASA  Headquarters 
Washington,  DC  205^6 

Attention:  C.  F.  Bersch,  RTM-6 

NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 

Attention:  Dr.  E.  Zaretsky  1 

W.  A.  Sanders  (49-1)  1 

and  Dr.  T.  Hergel  1 

Defense  Advanced  Research  Project  Office 
l400  Wilson  Blvd 
Arlington,  VA  22209 
Attention:  Dr.  Van  Reuth 

Inorganic  Materials  Division 
Institute  for  Materials  Research 
National  Bureau  of  Standard 
Washington,  DC  20234 

Department  of  Engineering 
University  of  California 
Los  Angeles,  CA  90024 

Attention:  Profs.  J.  W.  Knapp  and  G.  Sines 

Engineering  Experiment  Station 
Georgia  Institute  of  Technology 
Atlanta,  GA  30322 
Attexjtion:  J.  D.  Walton 

University  of  Illinois 
College  of  Engineering 
204  Ceramics  Building 
Urhana,  IL  6l801 
Attention:  Sherman  D.  Brown 

Department  of  Ceramic  Engineering 

Department  of  Engineering  Research 
North  Carolina  State  University 
Raleigh,  NC  27607 
Attention;  Dr.  H.  Palmour 

Ceramic  Science  and  Engineering  Section 
Pennsylvania  State  University 
University  Park,  PA  l6802 
Attention:  Dr.  R.  E.  Tressler 

Rensselaer  Polyteciiuic  Institute 
110  Eighth  Street 
Troy,  NY  12l8l 
Attention:  R.  J.  Diefendorf 


1 


1 


2 


1 


1 


1 


1 


1 


1 


1 
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NO  OF  COPIES 

School  of  Ceramics  1 

Rutgers,  The  State  University 
New  Brunswick,  NJ  08903 

Virginia  Polytechnic  Institute  1 

Minerals  Engineering 

Blacksburg,  VA  201+60 

Attention;  Dr.  D.  P.  H.  Hasselman 

Advanced  Mechanical  Teclmology  Inc.  1 

lUl  California  3:’fct 

Newton,  MA  02158 

Attention;  Dr.  Walter  D.  Syniuta 

Aerospace  Corporation  1 

Materials  Laboratory 

P.O.  Box  95085 

Los  Angeles,  CA  9001+5 

Batelle  Memorial  Institute  2 

505  King  Avenue 
Columbus,  OH  1+3201 

Attention;  Ceramics  Department  1 

Metal  &  Ceramic  Information  Center  1 

Research  and  Development  Division  1 

Carborundum  Company 
Niagra  Falls,  NY  11+302 
Attention;  Mr.  C.  McMurty 

Ceramic  Finishing  Company  1 

Box  1+98 

State  College,  PA  168OI 

Ceradyne  Inc .  1 

Box  1103 

Santa  Ana,  CA  92705 

Coors  Porcelain  Company  1 

600  Ninth  Street 

Golden,  CO  80I+OI 

Attention;  Research  Department 

Federal-Mogul  Corporation  1 

Anti-Friction  Bearing  R&D  Center 
3980  Research  Park  Drive 
Ann  Arbor,  MI  I+810I+ 

Attention;  D.  Glover 

Metallurgy  and  Ceramics  Research  Department  1 

General  Electric  R&D  Laboratories 
P.O.  Box  8 

Schenectady,  NY  12301 


NO  OF 

Hughes  Aircraft  Company 
Culver  City,  CA  90230 
Attention;  M.  N.  Gardos 

IIT  Research  Institute 
10  West  35th  Street 
Chicago,  IL  6o6l6 
Attention;  Ceramics  Division 

Kaweki-Berylco  Industry 
Box  lk62 

Reading,  PA  I9603 

Attention:  Mr.  R,  J.  Longnecker 

Research  &  Development  Division 
Arthur  D.  Little  Company 
Acorn  Park 
Cambridge,  MA  021  ^tO 

Mechanical  Technology,  Inc. 

968  Albany-Shaker  Road 
Latham,  NY  12110 
Attention:  Dr,  E.  F.  Finkin 

Norton  Company 
Industrial  Ceramics  Division 
One  New  Bond  Street 
Worchester,  MA  OI606 
Attention:  Dr.  M.  Torti 

Ceramic  Division 
Sandia  Corporation 
Albuquerque,  NM  87101 

SKF  Technology  Services 
SKF  Industries  Inc, 

1100  First  Street 
P.O.  Box  515 

King  Of  Prussia,  PA  19^06 
Attention:  P.  A.  Madden 

Southwest  Research  Institute 
P,  0.  Drawer  28 510 
San  Antonio,  TX  78228 

Materials  Sciences  &  Engineering  Laboratory 
Stanford  Research  Institute 
Menlo  Park,  CA  9^025 
Attention:  Dr.  Cubiciotti 
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Teledyne  CAE  1 

1330  Laskey  Road 
Toledo,  OH  U36OI 
Attention:  Hugh  Gaylord 

Union  Carbide  Corporation  1 

Parma  Technical  Center 
P.  0.  Box  6116 
Cleveland,  OH  UUlOl 

Materials  Sciences  Laboratory  1 

United  Technologies 
East  Hartford,  CT  061OI 
Attention:  Dr.  J.  J.  Brennan 

United  Technologies  Research  Center  1 

East  Hartford,  CT  06108 
Attention;  F.  L.  Ver  Snyder 

Astronuclear  Laboratory  1 

Westinghouse  Electric  Corporation 
Box  IO86I+ 

Pittsburgh,  PA  15236 

Westinghouse  Research  Laboratories  1 

1310  Beulah  Road 

Churchill  Borough 

Pittsburgh,  PA  15235 

Attention:  Dr.  R.  Bratton 


